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Image simulation techniques were employed to generate synthetic 
aperture radar (SAR) Images of a 17.7 km x 19.3 km test site located 
east of Lawrence, Kansas. The simulations were performed for a space 
SAR at an orbital altitude of 600 km, with the fol laving sensor para- 
meters; frequency ■ I»,75 GHz, polarization ■ HH, and angle of Inci- 
dence range « 7® - 22* from nadir. Three sets of Images were produced 
corresponding to three different spatial resolutions, namely 20 m x 
20 m with 12 looks, 100 m x 100 m with 23 looks, and 1 km x 1 km with 
1000 looks. Each set consisted of images for four different soil 
moisture distributions across the test site. The purpose of this 
study Is to evaluate the accuracy with which soil moisture can be pre- 
dicted for each of the 12 resolutlon/sol 1 moisture distribution combi- 
nations. The input Infonmatipn used to specify the gray level of each 
of the 800,000 pixels contained In the Image Included (when applicable) 
soil moisture, soil type, vegetation cover, surface roughness, row 
direction (relative to the radar look-direction), and local slope, 
while the prediction algorithm Is based on a generalized formula re- 
lating the received power to soil moisture, with no Information avail- 
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able to It on the scene properties except for the angle of Incidence 
with respect to the mean elevation of the test site. The results indicate 
that, for the agricultural portion of the test site, the soil moisture 
of about 90% of the pixels can be predicted with an accuracy of 120% 
of field capacity. Among the three spatial resolutions, the 1 km x 1 km 
resolution gave the best results for most cases; however, for very dry 
soli conditions, the 100 m x 100 m resolution was slightly superior. 


TABLE OF CONTENTS 


faae 

LIST OF FIGURES ...» ..... ill 

LIST OF TABLES xH 

ABSTRACT xv 

1.0 INTRODUCTION 1 

2.0 REVIEW OF THE RADAR RESPONSE TO SOIL MOISTURE 3 

2.1 Small-Scale Surface Roughness 5 

2.2 Periodic Surface Patterns 7 

2.3 Soil Texture . 21 

2.4 Vegetation Cover ............. 27 

2.5 Composite Effects of Scene Parameters . 31 

2.6 Summary of Sensor-Configuration Selection 36 

3.0 SYSTEM CONFIGURATION STUDY 39 

3.1 Introduction 39 

3.2 Equivalent Resolutions 43 

3.3 System Configuration Studies 45 

3.3.1 Real -Aperture Radar 45 

3.3.2 Synthetic Aperture Radar .. 50 

3.3.3 Scanning Radiometer-SAR 56 

3.4 System Tradeoffs 56 

3.5 Relation of Tradeoffs to Simulations 66 

4.0 SIMULATION STUDY 69 

4.1 Radar Image Simulation ..... 69 

4.1.1 An Overview of the Radar Image Formation 

Process 71 

4.1.2 Target Geometry Considerations .... 73 

4.1.3 Earth Curvature Considerations . 74 

4.1.4 Simulation of Fading and Other 

System Effects 75 

4.1.5 Noncoherent Processing to Increase Averaging . . 80 

4.1.6 Radar System Parameters Modeled 80 

4.2 Determination of Mean Backscattering Coefficient of 

Target Conditions as a Function of Incidence 

Angle and Moisture ..... 81 

i 



Paje 


4.2.1 Bare-Soil Algorithms , ...... 

4.2.2 Vegetation-Covered Soil Algorithms ....... 

4.2.3 Target Classes with No Dependence on 

Soil Moisture 

4.2.4 Angular Dependence of 5° for Given 

Moisture Conditions . . 

4.3 Simulation Data-Base Construction ........... 

4.3.1 Target Categories <• . 

4.3.2 Soil Textural Classification . 

4.3.3 Surface Elevation 

4.3.4 Data Base Registration ....... 

4.4 Simulation of Slope , • 

4.5 Simulation of Rainfall and Soil Moisture Conditions . . 

4.5.1 Soil Moisture Conditions 

4.5.2 Postulated Rainfall and Evaporative History . . 

4.6 Generation of Radar Image— Summary of Computer 

Algorithms and Images . 

5.0 EVALUATION OF RESULTS . 

5.1 Generalized Moisture Interpretation Algorithms . . . . 

5.2 Interpretation of Simulated Radar Images ....... 

5.3 Evaluation of Moisture Estimate Accuracy Using 

Blind Classifier . 

5.4 Analysis of Moisture Estimate Accuracy Over 

the Total Data Base 

5.5 Analysis of Moisture Estimate Accuracy 

Within the Floodplain . 

Effect^! of Crop Type, Surface Roughness 
and Row Direction 

6.0 CONCLUDING REMARKS 

REFERENCES 


87 

90 

92 

93 

93 

96 

101 

104 

106 

106 

108 

108 

no 

115 

118 

126 

140 

142 

151 

166 

175 

182 

185 


ii 


LIST OF FIGURES 


Figure 2.1 


Figure 2.2 


Figure 2.3 


Figure 2.4 
Figure 2.5 


Figure 2.6 


Figure 2.7 


Figure 2.8 


Figure 2.9 


Page 


Linear correlation coefficient between 
o°(dB) and volumetric soil moisture of 
the surface 0-1 cm layer, plotted as a 
function of angle of Incidence, for five 
different surface roughnesses (from [4]) ... 6 

Linear correlation coefficient between 

o°(dB) and volumetric soil moisture as 

a function of angle of Incidence for 

L-band (1.6 GHz), C-band (4.75 GHz), 

and X-band (0.4 GHz) (from Jackson 

et al. [17]) 8 

Variation of a° with RMS height at 4.5 GHz, 

HH polarization for several angles of 
Incidence. Soil moisture range: 0.12 to 
0,18 g/cm® (from Le loan et al. [8]) 9 

Variation of a® with volumetric soil moisture 
at (a) 1.5 GHz, (b) 3 GHz and (c) 4.5 GHz 
(from Le loan et al. [8]) 10 


Comparison of the angular response of the 

look-direction modulation function of a 

com field for HH, HV, and VV polarizations 

at (a) 1.1 GHz and (b) 4.25 GHz 

(from [13]) 15 

Comparison of the angular response of the 
look-direction modulation function of a 
soybean field for HH, HV, and VV polariza- 
tions at (a) 1.1 GHz and (b) 4.25 GHz 
(from [13]) 16 

Cwiparison of the angular response of the 

look-direction modulation function of a 

wheat field for HH, HV, and VV polarizations 

at (a) 1.1 GHz and (b) 4.25 GHz 

(from [13]) .17 


Look-direction modulation function as a 
function of angle of incidence for a row- 
spacing of 100 cm and row height of 25.4 cm 


(from Fenner et al. [10]) ..... 18 

Look-direction modulation function as a 
function of angle of Incidence for a row- 
spacing of 100 cm and row height of 15 cm 
(from Fenner et al. [10]) . .19 


11i 


Figure 2.10 

Figure 2.11 
Figure 2.12 
Figure 2.13 
Figure 2.14 

Figure 2.15 
Figure 2.16 
Figure 2.17 


Pace 

tmmrnimmm 


Variation of o° with angle of Incidence for a 
wheat- stubble field, from airborne observation 
flown In parallel and perpendicular directions 
relative to row direction. It Is suspected that 
the observed differences between a // ando°i for 
HV polarization are not real; theaoJ_data are 
biased due to polarization coupling by the antenna 
(from Cl8] ) . 20 


Dielectric constant as a function of volumetric 
moisture content at 5 GHz for the soils listed 
In Table 2.2 23 

Dielectric constant as a function of volumetric 

moisture content at 1,412 GHz for the soils 

listed In Table 2.2 (fromC22]) . 24 

Reflection coefficient R(dB) of sand, clay loam, 
and clay at 1.4 GHz, 0 degrees as a function of 
volumetric moisture (from [5)) 25 


Backscattering coefficient at 4.6 GHz, 10 degrees, 

HH polarization as a function of 0-1 cm gravimetric 
soli moisture (a) sandy loam, (b) silty clay loam, 
and (c) silty clay. Saturated soil conditions In 
the field of silty clay loam are not Included In 
linear regression (from [5]) 26 

Linear regression fits of radar response at 

4.625 GHz, 10 degrees, and HH polarization to 

volumetric soil moisture In the 0-1, 0-2, 0-5, 

and 0-9 cm layers. Regression results are 

based on all data obtained for each soil texture 

and an equally distributed combined texture 

data base (from [5]) 28 


Linear regression fits of radar response at 
4.625 GHz, 10 degrees, and HH polarization 
to percent of 0.33 bar water content In the 
• 0-1, 0-2, 0-5, and 0-9 cm soil layers. 

Regression results are based on all data obtained 
for each soil texture and an equally distributed 
combined texture data base (from [5]) 29 

Linear correlation coefficient as a function of 
frequency for var^lous 0-5 cm soil moisture 
Indicators with a at 10 degrees incidence angle 
and HH polarization. Soil moisture is expressed 
as gravimetric Mg, volumetric My, percent of 0.33 
bar moisture, percent of estimated field capacity 
Nfc, and percent of 1.0 bar moisture. Data 
1s from the combined multi texture data base 
(from [5]) . . ^ 30 


1v 


Figure 2.18 
Figure 2.19 
Figure 2.20 
Figure 2.21 
Figure 2.22 

Figur; 2.23 


Figure 3.1 

Figure 3.2 
Figure 3.3 
Figure 3.4 
Figure 3.5 


Figure 4.1 


Soil moisture content In percent of field 
capacity as a function of a®(dB) at 4.25 GHz, 

0 ■ 10® and HH polarization (from [14]) ..... 32 


Radar backscatter response to soil moisture 
for bare fields, vegetation-covered fields, 
and both types combined 33 

Radar backscatter response to soil moisture 
based on airborne measurements at 4.75 GHz, 

0 ■ 10® and HH polarization (from [18]) ..... 34 


Variation of a° with percent field capacity 
for different covers, at 4.75 GHz, HH polari- 
zation and e « 1^)®, (from [18]) 35 

Correlation coefficient, sensitivity and 
Y-Intercept for linear regression of a®(dB) 
against percent field capacity, at 4.75 GHz, 


HH polarization (from [18]) ........... 37 

Correlation coefficient, sensitivity and 
Y-intercept for linear regression of a°{dB) 
against percent field capacity, at 4.75 GHz, 

HH polarization (from [18]) % . . . 38 

Side-looking SAR observing a curved earth 

from an altitude of 600 km 42 

Configuration of RAR . 46 

SAR configuration assumed ... 51 

RADISAR configuration 57 


(a) Average transmitter power, and (b) telemetry 
rate, as a function of equivalent-square 
photographic-quality pixel dimension. The numeral 
next to a point refers to the case number 
In Table 3.1 64 

Backscattering coefficient a° at 4.5 GHz, 

HH polarization, 10® angle of incidence from smooth 
bare soil plotted as a function of percent of 
1/3-bar water content In the 0-5 cm 
soil layer 84 


V 


Figure 4.2 Mean backscattering coefficient of several 
target classes as a function of angle of 
Incidence for very dry soil moisture 
conditions, H|;q ^ 10,0!K ............ 94 

Figure 4.3 Mean backscattering coefficient of several 
target classes as a function of angle of 
Incidence for saturated soil moisture 
conditions, MpQ ** ?50.0%. 95 

Figure 4.4 High altitude Infrared Imagery obtained 
May 15, 1978 and used to classify land** 
use categories and field boundaries within 
the data base . ....... 97 


Figure 4.5 Location of non-agri cultural features 

within the simulation data base. Trees, 
highways, and secondary roads are shown 
In black; water bodies, railroads, bridges, 


and buildings are shown In white. ....... 97 

Figure 4.6 Location of bare soil fields within the 
simulation data base. Random surface 
roughness conditions from smooth to rough 
are proportional to image intensity ...... 98 

Figure 4.7 Location of vegetated fields within the 
simulation data base. Hay, pasture, 
soybeans, milo, corn, and wheat are shown 
with ascending image intensity 98 


Figure 4.8 Soil classes within the simulation 

test site. Pixel intensity is related 
to a horizon soil texture from 
Table 4.6 where sand is darkest and 
the complexes are brightest. The 
Image represents a sampling of every 
other column and ravi of the entire data 


base (25!^ sample) 103 

Figure 4.9 Elevation above mean sea level. Every 
other column and row of the total data 
matrix are sampled to produce this 
Image 105 

Figure 4.10 Input soil moisture for moisture 
condition 1, three days after a 
saturating rain. All soils are at 
100% of field capacity. . * 113 


v1 


Figure 4.11 

Figure 4.12 
Figure 4.13 

Figure 4.14 

Figure 4,15 

Figure 4.16 

Figure 4.17 

Figure 4.18 

Figure 4.19 

Figure 4.20 

Figure 4.21 


Input soil moisture for Moisture Condition 2, 
Immediately after passage of a thunderstorm* 


Surface soil moisture ranges from field 

capacity to saturation ...*,* 113 

Input soil moisture for Moisture Condition 3, 

10 days after the thunderstorm ... 114 

Input soil moisture for Moisture Condition 4, 

30 days after the thunderstorm, 

.drought conditions . 114 


Radar simulation with 20-meter resolution 

and 7.5° angle of Incidence for Moisture 

Condition 1, all soils a't field capacity .... 120 


Radar simulation with 20-meter resolution 
and 7.5° angle of incidence for Moisture 
Condition 2, Immediately after 
a thunderstorm ........ 120 

Radar simulation with 20-meter resolution 

and 7.5° angle of Incidence for 

Moisture Condition 3, 10 days after 

the thunderstorm .... 121 

<» 

Radar simulation with 20-meter resolution 

and 7.5° angle of Incidence for 

Moisture Condition 4, 30 days after 

the thunderstorm .... 121 

Radar simulation with 100-meter resolution 

and 7.5° angle of Incidence for 

Moisture Condition 3. 10 days after 

the thunderstorm 122 

Radar simulation with l-k11ometer resolution 

and 7.5° angle of Incidence for 

Moisture Condition 3, 10 days after 

the thunderstorm .... ... 122 

Radar simulation with 20-meter resolution 

and 11.1° angle of Incidence for 

Moisture Condition 3, 10 days after 

the thunderstorm 123 

Radar simulation with 100-meter resolution 

and 11.1° angle of Incidence for 

Moisture Condition 3, 10 days after 

the thunderstorm 123 


v11 


Figure 4.22 


Figure 5.1 


Figure 5.2 


Figure 5.3 


Figure 5.4 


Figure 5.5 


Figure 5.6 


Figure 5.7 


Figure 5.8 


Figure 5.9 


Rader slimiUtlon with l-k11(xneter 

resolution and 11. 1** angle of Incidence 

for Moisture Condition 3, 10 days after 

the thunderstorm ........ 124 

Backicatter coefficient as a function of Mpc 
In the 0-5 on soil layer for 324 observations 
of bare and vegetated soils at 4.5 GHz, 

HH polarization, and 10** Incidence angle .... 129 

Backscattering coefficient as a function of Mpc 
In the 0-5 cm soil layer for 324 observations 
of bare and vegetated soils at 4.5 GHz, 

HH polarization and 20® angle of Incidence ... 130 

Backscattering coefficient as a function of Mpc 
In the 0-5 cm soil layer for 181 observations 
of bare soil at 4.5 GHz, HH polarization and 
10® angle of Incidence . 131 

Backscattering coefficient as a function of Mpc 
In the 0-5 cm soil layer for 181 observations 
of hare soil at 4.5 GHz, HH polarization, and 
20® angle of incidence ............. 132 

Backscattering coefficient as a function of Mpc 
In the 0-5 cm layer for 143 observations of 
various crops at 4.5 GHz, HH polarization, and 
10® angle of Incidence 134 

Backscattering coefficient as a function of Mpc 
In the 0-5 cm layer for 143 observations of 
various crops at 4.5 GHz, HH polarization, and 
20® angle of Incidence 135 

Intercept of generalized algorithms as a 
function of angle of incidence between 0® and 
30®. f(0) Is shown for (a) all agricultural- 
scene data combined, (b) bare soil data only, 
and (c) crop-canopy data only 138 


Slope of generalized algorithms as a function 
of angle of Incidence between 0® and 30®. 
g(e) is shown for (a) all agricultural^scene 
data combined, (b) bare soil data only, and 
(c) crop-canopy data only ........... 139 

General Interpretation procedure for evaluating 
soil moisture estimate accuracy 141 


vlii 


Figure 5.10 


Figure 5.11 


Figure 5.12 


Figure 5.13 


Figure 5.14 


Figure 5,15 


Figure 5.16 


Figure 5.17 


Figure 5.18 


Figure 5.19 


Predicted moisture rwp estimated by 
general altjoHthm from radar simulation 
with 20-meter resolution and 7,5“ angle 
of Incidence for Moisture Condition 1 , , , . 143 

Predicted moisture map estimated by 
general algorithm from radar simulation 
with 20-meter resolution and 7,5® angle 
of Incidence for Moisture Condition 2 , , ♦ . 143 

Predicted moisture map estimated by 
general algorithm from radar simulation 
with 20-meter resolution and 7.5° angle 
of Incidence for Moisture Condition 3 . . . . 144 

Predicted moisture map estimated by 
general algorithm from radar simulation 
with 20-meter resolution and 7,5° angle 
of Incidence for Moisture Condition 4 , . . . 144 

Predicted moisture map estimated by 
general algorithm from radar simulation 
with 100-meter resolution and 7,5° angle 
of Incidence for Moisture Condition 3 . . , . 145 

Predicted moisture map estimated by 
general algorithm from radar simulation 
with 1-kilometer resolutii,i, and 7,5° angle 
of Incidence for Moisture Condition 3 , , . . 145 

Soil moisture estimate error, predicted Mpc 
minus Input Mfc. from radar simulation with 
20-meter resolution and 7,5° angle of 
Incidence for Moisture Condition 3 148 

Soil moisture estimate error, predicted Mpp 
minus Input Mpc, from radar simulation with 
100-meter resolution and 7.5° angle of 
Incidence for Moisture Condition 3 . . , , 148 

Soil moisture estimate error, predicted Hpc 
minus Input Mpc* from radar simulation with 
l-kllometer resolution and 7.5° angle of 
Incidence for Moisture Condition 3 .... 149 

Soil moisture estimate error, predicted Mpc 
minus input Mpc, from radar simulation 
with 100-meter resolution and 11.1° angle 
of Incidence for Moisture Condition 3 . . . . 149 


ix 


Figure 5.20 


Figure 5.21 


Figure 5.22 


Figure 5.23 


Figure 5.24 


Figure 5.25 


Figure 5.26 


Figure 5.27 


Figure 5.28 


Figure 5.29 


Soil liroisture estimate error level -si iced 

at Accuracy Levels 3 and 4 from radar 

simulation with 20-Fneter resolution and 

7.5® angle of incidence for Moisture 

Condition 3 152 

Soil moisture estimate error level>-sliced 

at Accuracy Levels 3 and 4 from radar 

simulation with 100-meter resolution and 

7.5® angle of incidence for Moisture 

Condition 3 ...... 152 

Soil moisture ertimate error level -si iced 

at Accuracy Levels 3 and 4 from radar 

simulation with 1-kilometer resolution and 

7.5® angle of incidence for Moisture 

Condition 3 .153 

Soil moisture estimate accuracy for the 

entire simulation data base with Moisture 

Condition 1, all soil at field capacity, and 

at angle of incidence between 7.5® and 9.3® .... .159 

Soil moisture estimate accuracy for the entire 
simulation data base with Moisture Condition 2, 
saturated soil, and at angles of incidence between 


7.5® and 9.3® 160 

Soil moisture-estimate accuracy for the entire 
simulation data base with Moisture Condition 3, 

10 days after thunderstorm, and at angles of 
incidence between 7.5® and 9.3® 161 

Soil moisture estimate accuracy for the entire 
simulation data base with Moisture Condition 4, 
drought, and at angles of incidence between 
7.5® and 9.3® .162 

Soil moisture estimate accuracy for the entire 
simulation data base with Moisture Condition 3, 

10 days after thunderstorm, and at angles of 
incidence between 11.1® and 12.8® .... .163 


Soil moisture estimate accuracy for agricultural 
categories on the floodplain with Moisture 
Condition 1, all soil at field capacity, and 
at angles of incidence between 7.5® and 9.3° 170 

Foil moisture estimate accuracy for agricultural 
categories on the floodplain with Moisture 
Condition 2, saturated soil, and at angles of 
incidence between 7.5® and 9.3® .171 


X 



Figure 5.30 


Figure 5.31 


Figure 5.32 


Figure 5.33 


Figure 5,34 


Soil moisture estimate accuracy for 

agricultural categories on the floodplain 

with Moisture Condition 3, 10 days after 

a thunderstorm, and at angles of Incidence 

between 7,5° and 9,3° . 172 

Soil moisture estimate accuracy for 

agricultural categories on the floodplain 

with Moisture Condition 4, drought, and at 

angles of incidence between 7,5° and 9.3° . * . . . 173 

Soil moisture estimate accuracy for 

agricultural categories on the floodplain 

with Moisture Condition 3, 10 days after a 

thunderstorm and with angles of incidence 

between 11.1° and 12.8° 174 

Percent of data base where estimate moisture 
is within ±30^ of actual moisture Mpr 
Angle of incidence 6 is 7,5° to 9.3° relative 
to nadir. Accuracy level 3 is approximately 
±30% of field capacity, . 176 

Percent of data base where estimated moisture ^^c 
Is within ±20% of actual moisture Mpg. Angle 
of Incidence 9 Is 7.5° to 9.3° relative to nadir. 
Accuracy level 4 is approximately ±20° of field 
capacity. . 177 


LIST OF TABLES 


Page 


TABLE 2.1 

Summary of Active Microwave Investigations 
of Soil Moisture and Related 
Scene Variables 

. 4 

TABLE 2.2 

Different Types of Soils Used for the Dielectric 
Measurements as a Function of Water Content 
at 1.4 GHz and 5 GHz (from [22]). 

. 22 

TABLE 3.1 

Design Options— Condensed Table . 

. 59 

TABLE 3.2 

Potential System Design Options 

. 62 

TABLE 3.3 

Data-Handling Parameters. 

. 63 

TABLE 3.4 

Simulations Compared with Nearest- 
System Case . . 

. 68 

TABLE 4.1 

Linear Regression Coefficients, Correlation 
Coefficients, Sample Size, and Data Sources 
of 4.75 GHz, HH Radar Response to 0-5 cm 
Soil Moisture of Target Glasses at 
Selected Incidence Angles ...... 

. 83 

TABLE 4.2 

Mean Algorithms for Simulated 

Target Classes. 

. 85 

TABLE 4.3 

RMS Surface Roughness and Soil Textural 
Classification of Experimental Data 
Used to Develop Bare-Soil Algorithms. ..... 

. 89 

TABLE 4.4 

Area Percent of Total Data Base 
Assigned to Each Target Class . 

. 99 

TABLE 4.5 

Comparison of Data Base Composition with 
Historical Records for Douglas 
County, Kansas 

.100 

TABLE 4.6 

Soil Textural Classes, Their Area Percent, 
and Characteristic 1 /3-Bar Water Contents 
Within the Simulation Data Base ... 

.102 

TABLE 4.7 

Values of Rainfall and Evaporation Used 
to Derive the Hypothetical Soil 
Moisture Conditions ....... 

.112 

TABLE 4.8 

Summary of Simulated Radar Images . . 

.119 

TABLE 5.1 

Results of Least-Square Linear Regressions 
Used to Establish the General Interpretation 
Algorithms for (1) All Classes, (2) Bare 
Soil Classes, and (3) Crop-Canopy Classes . . . 

.128 



xii 



Page 


TABLE 5.2 Third-Order Polynomial Fits of Linear 
Regression Coefficients as a Function 
of Incidence Angle Between 0° and 30° 
for (a) All Agricultural Scene Data, 
(b) Bare Soil Data Only, and (c) Crop 
Canopy Data Only 


TABLE 5.3 Interpretation of Soil Moisture Estimate 

Error Maps . . . * 154 

TABLE 5.4 Soil Moisture Information and Data 
Requirements at Different Crop 
Production Stages (from [15]) 155 

TABLE 5.5 Needs of Soil Moisture Information and 
Data Requirements in Hydrology 

(from [51]) 156 

TABLE 5.6 Approximate Relationship Between 
Uncertainty in Percent of Field 
Capacity and User Requirements of 
Moisture Estimate Accuracy 157 

TABLE 5.7 Cumulative Percent of Total Data Base 

Within Certain Absolute Estimate Error 

Limits for C-Band Simulation Results 

at e = 7.5° - 9.3°. 164 

TABLE 5.8 Cumulative Percent of Total Data Base 

Within Certain Absolute Estimate Error 

Limits for C-Band Simulation Results 

at 0 = 11.1° - 12.8° 165 

TABLE 5.9 Cumulative Percent of Agricultural Flood- 
plain Within Certain Absolute Error 
Limits for C-Band Simulation Results 
at 9 = 7.5° - 9.3° 168 

TABLE 5.10 Cumulative Percent of Agricultural 
Floodplain Within Certain Absolute 
Error Limits for C-Band Simulation 
Results at e * 11.1° - 12.8° 169 

TABLE 5.11 Area Percent of the Agricultural Flood- 
plain for Crop Categories. Floodplain 
is Defined by Elevation $ 820 Feet 178 

TABLE 5.12 ±2035 Accuracy Levels for Crop Categories 

on the Floodplain for the 7.5“ to 9.3° 

Angular Swath 179 


xiii 




TABLE 5.13 i20X Accuracy Levels for Crop Categories 
on the Floodplain for the 11.1® to 12.8® 
Angular Swath .... 



EVALUATION OF THE SOIL MOISTURE PREDICTION 
ACCURACY OF A SPACE RADAR USING 
SIMULATION TECHNIQUES 


F. T. Ulaby, C. Dob$on» 0. Stiles, R. K, Moore and J. Holtztnan 

Remote Sensing Laboratory 
University of Kansas Center for Research* Inc. 

Lawrence, Kansas 66045 

ABSTRACT 

Image simulation techniques were employed to generate synthetic 
aperture radar (SAR) images of a 17.7 km x 19.3 km test site located 
east of Lawrence, Kansas. The simulations were performed for a space 
SAR at an orbital altitude of 600 km, with the following sensor para- 
meters: frequency * 4.75 GHz, polarization « HH, and angle of inci- 
dence range = 7° - 22° from nadir. Three sets of Images were produced 
corresponding to three different spatial resolutions, namely 20 m x 20 m 
with 12 looks, 100 in x 100 m with 23 looks, and 1 km x 1 km with 1000 
looks. Each set consisted of images for four different soil moisture 
distributions across the test site, The purpose of this study is to 
evaluate the accuracy with which soil moisture can be predicted for each 
of the 12 resol ution/soll moisture distribution combinations. The Input 
Information used to specify the gray level of each of the 800,000 pixels 
contained in the image included (when applicable) soil moisture, soil 
type, vegetation cover, surface roughness, row direction (relative to 
the radar look direction), and local slope, while the prediction algo- 
rithm is based on a generalized formula relating the received power to 
soil moisture, with no Information available to It on the scene pro- 
perties except for the angle of incidence with respect to the mean 
elevation of the test site. The results indicate that, for the agri- 
cultural portion of the test site, the soil moisture of about 9055 of 
the pixels can be predicted with an accuracy of ±2056 of field capacity. 
Among the three spatial resolutions, the 1 km x 1 km resolution gave 
the best results for most cases; however, for very dry soil conditions, 
the 100 m x 100 m resolution was slightly superior. 
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1.0 INTRODUCTION 


During the past decade, several Investigations have been 
conducted to evaluate the role of active microwave remote sensing 
techniques for sensing soil moisture content [1-19]. Although 
some of these investigations were not strictly related to soil 
moisture, their findings provided Information on scene features 
(such as vegetation cover and surface roughness) that are related 
to the problem of sensing soil moisture with radar. 

With the exception of one study, the Investigations reported 
to date have Involved the use of ground-based [1-14] or airborne 
radars [10, 13-18]. The exception Is a Skylab Investigation [19] 

In which the backscattering coefficient measured by Skylab 's 
13.9 GHz scatterometer, was related to the moisture content of 
soil samples obtained from the ground areas observed by the sensor. 
Although a° exhibf*^;ed a fair degree of correlation with moisture 
content, the results of the Investigation should be considered as 
only a positive Indication, rather than as a demonstration, of the 
potential use of radar for sensing soil moisture. This qualification 
Is based on the fact that the resolution area of Skylab 's scatterometer 
was of the order of 200 km*, while the soil moisture estimate was based 
on a few point-samples. 

Seasat provided the first opportunity for Investigating the 
radar response to soil moisture with a high-resolution spaceborne 
Imaging radar. In support of a project organized by Texas A&M 
University, Seasat SAR imagery was acquired for several passes over 
a soil moisture test-site in Oklahoma, and aircraft underflights were 
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also made using the NASA/ JSC 1.6 GHz scatterometer, A report sum- 
marizing the results Is forthcoming [20]. 

The purpose of this report is to document the results of an 
investigation conducted tp evaluate the applicability of active 
microwave approaches to soil -moisture sensing from a satellite platform 
and to define the general system parameters of such an active microwave 
system. The report is divided into three major topics. The first 
topic (Section 2) consists of a literature review of the radar 
response to soil moisture and related scene parameters. The objective 
of this section is to specify the radar configuration that holds the 
greatest promise in terms of the accuracy and precision with which soil 
moisture may be estimated with an active microwave system. Specifically, 
the following parameters are to be identified: (a) microwave frequency, 

(b) angle of incidence range, (c) polarization configuration, and 
(d) the lower and upper limits of the dynamic range of o°. These para- 
meters are used in Section 3 to evaluate the impact of spatial resolution 
on system-hardware and signal -processing requirements. Several different 
system configurations are considered, including the fully focused 

synthetic aperture radar (SAR), partially focused SAR, unfocused SAR, 

* 

and the real aperture radar (RAR). Additionally, a combined radiometer- 
SAR (RADISAR) configuration is considered. 

Three candidate configurations will be chosen to evaluate the 
relationship between soil moisture estimation accuracy and spatial 
resolution. This evaluation is the subject of the third topic 
(Sections 4 and 5). For each candidate configuration, image simula- 
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tion teahniques are employed to generate a radar Image of a test site 
located In the vicinity of Lawrence, Kansas, using a multitude of 
scene parameters to characterize the scattering coefficient of each 
resolution cell (pixel). Among these parameters are the soil mois- 
ture content, soil texture, local slope, cover-type, etc. The image 

thus generated is then subjected to a generalized algorithm using 
0 

a and satellite altitude and position as sole input^^ >?nd moisture 
content as sole output. Upon applying this algorithm, a predicted 
soil moisture "map" is produced. A soil moisture error map is then 
produced by taking the difference between the predicted moisture 
content and the true moisture content on a pixel -by-pixel basis. 

The performance of each of the three candidate system configurations 
(corresponding to three different spatial resolutions) is evaluated 
in terms of error distribution functions. 

2.0 REVIEW OF THE RADAR RESPONSE TO SOIL MOISTURE 

Table 2.1 provides a summary of the major studies reported in 
the literature that pertain to radar backscatter from soil surfaces. 
The scene parameters of primary importance are: 

(a) soil moisture (profile) 

(b) small-scale soil surface roughness 

(c) large-scale periodic surface patterns (row tillage) 



Sumaiy of Active Hfcrtwave Investigations 
of Soil Moisture and Related Scene Variables 
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(d) soil texture 

(e) vegetation cover 

Parameter (a), soil rr,3lsture content, Is the quantity a interest, 
while parameters (b) • (e) are quantities whose Influence on the 
backscattering coefficient should be either minimized or accounted 
for, as discussed next. 

2.1 Small-Scale Surface Roughness 

Surface roughness Is the prime mechanism responsible for making 
non-nadir radar observations possible, since for a perfectly smooth 
surface the backscattered pow'jr Is zero except at nadir. Hence, the 
effect of surface roughness on 0 *^ has received the greatest attention, 
relative to the other factors listed above. Based on experiments 
conducted by the University of Kansas in 1974 (for three different 
surface roughnesses) and in 1975 (for five different surface rough- 
nesses), a set of radar parameters was proposed for minimizing the 
effect of surface roughness on 0 °. These were: 4GHz 5 : f 4 : 5GHz, 

7® < 0 < 20® and HH -polarization. Comparable performance may also 
be achieved using HV polarization, but the use of HV polarization Is 
usually less desirable than HH because of the higher transmitter 
power levels that are required. Figure 2.1 shows a summary of the 
results In the form of the correlation coefficient p as a function 
of 0 for three different microwave frequencies. Here p Is the linear 
correlation coefficient between a°(dB) and the volumetric moisture 
content of the surface 0-1 cm layer. Later analyses have shown that 
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Figure 2.1 Linear correlation coefficient between o°(dB) and volumetric soil 

moisture of the surface 0-1 cm layer, plotted as a function of angle 
of incidence, for five different surface roughnesses (from [4]). 
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slightly higher magnitudes of p are obtained If the correlation Is 
performed against the moisture In the 0-5 cm soil layer. 

Figure 2.2 shows a similar set of curves for p versus e based on 
airborne measurements made over watersheds In Oklahoma [17]. Although 

h 

the magnitude of p Is generally lower than that shown In Figure 2.1, 
the overall angular and frequency behavior Is very similar. The lower 
magnitude of p Is attributed to surface slope effects. The C-band 
channel (at 4.75 GHz) provides the highest magnitude of p (among the 
three frequencies used) and the peak value of p Is In the 10® - 15® 
region, In agreement with the results shown In Figure 2.1. 

The above conclusions were confirmed further by the experimental 
results reported by Le Toan et al. [8] of the Paul Sabatier University 
In Toulouse, France. Samples of their results are shown In Figures 
2.3 and 2.4. The dependence on surface roughness Is Illustrated In 
Figure 2.3 where o°(dB) at 4.5 GHz Is plotted as a function of rm^ 
surface height for several angles of Incidence e. On the basis of the 
behavior depicted In Figure 2.3 and other analyses, it was concluded 
that a° Is least dependent on rms height for e in the 6° - 20® range. 
Figure 2.4 shows versus m^, the volumetric moisture content, at 
three microwave frequencies. The data shown Include the smoothest and 
roughest surface conditions observed in their investigation. 

2 . 2 Periodic Surface Patterns 

Row tillage patterns associated with row crops, such as wheat, 
corn, and soybeans, represent a large-scale periodic surface pattern, 
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Figure 2.4 
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Variation of a° with wlumetric soil moisture at (a) 1.5 GHz, 
(b) 3 GHz and (c) 4.5 GHz (from Le loan et al. [8]). 
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superimposed on which is the small-scale roughness discussed in the 
previous section. For rectangular fields with parallel rows, the peri- 
odic pattern 'iS unidirectional exhibiting a repetitive spatial modu- 
lation function in the direction perpendicular to the row direction. 

For fields planted in circular rows, the spatial modulation function is 
a function of r, the distance from the center of the field. 

The effect of periodic patterns on a° has be^* examined by three 
experimental investigations [10, 13, 18], These investigations were 
performed for two classes of periodic patterns; the observations re- 
ported by Ulaby and Bare [13] are for row spacings and depths that 

are typical of dry-land farming practices, whereas the other two in- 
vestigations [10, 18] are for comparatively deeper rows that are 
typical of irrigated-field conditions. Since the present study is 
based on conditions for a test site in eastern Kansas, where irrigation 
is seldom used, the results obtained by Ulaby and Bare [13] will be 
examined first. 

Although the Ulaby and Bare [IS- study was conducted for fields 
planted in wheat, soybeans, and corn, experimental evidence has shown 
that at frequencies below 8 GHz, the vegetation cover exercises a 
minor effect on at angles in the 0° - 20® range. Two types of 
periodic patterns were observed. The corn and soybean fields were 
characterized by a row spacing of 90 cm and a row height (trough to 
peak) of 6 cm, while the wheat fields had 15-cm row spacing and 2-cm 
row height. The difference between and a®^(d.B), corresponding 
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to observations with the look direction orthogonal and parallel to the 
row direction, respectively, has been named the look-direction modula- 
tion function, M(dB). Figures 2.5 - 2.7 show M(dB) as a function of 
«ngle for HH, HV and VV polarization. For the two frequencies 
shown, It is observed that M(dB) Is much higher at 1.1 GHz In compari- 
son to 4.25 GHz. Furthermore, generally HV polarization appears to be 
considerably less sensitive to look direction than the like-polarized 
configurations. On the basis of this study, it was concluded that for 
soil moisture determinations, the look-direction ambiguity problem may 
be avoided by operating at frequencies higher than 4 GHz with HH polar- 
ization, or at any frequency with HV polarization. It should be em- 
phasized that these conclusions pertain to the periodicities and row 
depths investigated in the above study, which are typical of dry-land 
farming soil surface conditions. 

For periodic patterns with deeper rows, a° was observed to be 
sensitive to look direction (relative to row direction) at 1.6 GHz, 

4.75 GHz and 13.3 GHz [10, 18]. Figures 2.8 and 2.9 show plots of 
M(dB) as a function of e for HH and VV polarizations, based on air- 
borne scatterometer data acquired by NASA/ JSC [10]. It is observed 
that the location of the peak is Independent of the frequency and that the 
magnitude decreases with increasing frequency. These results are in 
agreement with theoretical model predictions [21]. Another set of 
row-direction data Is shown In Figure 2.10 for a wheat stubble 
field [18]. Large differences are observed between and for 
HH polarization at all the frequencies. In contrast, smaller differences 
are observed for the HV-polarIzed data. Further discussion of this 
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Figure 2.5 Comparison of the angular response of the look- direction modulation 
function of a corn field for HV and VV polarizations at a) 1.1 
GHz and b) 4.25 GHz (from [13]). 


' Crop Type: Soybeans 
Crop Height (cm): 70 
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Figure 2.7 Comparison of the angular response of the look- direction modulation 
function of a wheat field for HH, HV and VV polarizations at a) 

1.1 GHz and b) 4.25 GHz {from [13]). 
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Figure 2.8 Look ‘Direction Modulation Function as a function of angle 
of incidence for a row spacing of 100 cm and row height 
of 25.4 cm (from Fenner et al. [10]). 
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of incidence for a row spacing of 100 cm and row height 
of 15 cm (from Fenner et al. [10]). 
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Figure 2.10 Variation of a° with angle of incidence for a wheat-stubble field, I 

from airborne observation flown in parallel and perpendicular directions | 

relative to row direction. It is, suspected that the observed differences | 

between o// and for HV polarization are not real; the o^ data are biased | 
due to polarization coupling by the antenna (from [18]). I 


subject Is deferred to Section 2.5 where the composite Influence of all 
the scene variables Is considered, 

2.3 Soil Texture 

Soil texture exercises an Indirect, but significant, Influence on 
the sensitivity of a° to soil moisture. Experimental measurements In- 
dicate that the dielectric properties of soil are dependent on Its 
composition. For a given volumetric moisture content, the real part 
of the dielectric constant is generally highest for sand and lowest 
for clay. Figures 2.11 and 2.12 contain a summary ot measured dielec- 
tric-constant curves for a variety of soils (Table 2,2) at 5 GHz and 

1.4 GHz, Generally, the spread In values due to texture decreases 
slowly with frequency. 

The scattering coefficient a° Is related to the reflection co- 
efficient, which in turn is governed by the dielectric properties of 
the soil. Thus, the reflection coefficient of different soil types 
exhibit different responses to soil moisture content. This is illus- 
trated in Figure 2.13 for 1.4 GHz (based on data reported by Newton 
[24]). The dependence of a° on soil texture Is shown In Figure 2.14 
[5]. 

The dependence on soil type may be removed by expressing soil 
moisture in terms of a tension-related quantity. Schmugge [25] has 
proposed the use of percent of field capacity, MpQ, as a soil moisture 
Indicator where field capacity is the gravimetric moisture content 
at 1/3-bar tension. Additionally, Schmugge provides an approximate 
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Table 2.2 Different types of soils used for the dielectric measurements as a 
function of water content at 1.4 GHz and 5 GHz (from [221. 
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Figure 2.11 Dielectric constant as a function of volumetric moisture 

content at 5 GHz for the soils listed in Table 2.2 (from [22]). 
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Figure 2,12 Dielectric constant as a function of volumetric moisture 
content at 1.412 GHz for the soils listed in Table 2.2 
(from [22]). 
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Figure 2.13 Reflection coefficient R(dB) of sand, clay loam, and clay 
at 1.4 GHz, 0 degrees as a function of volumetric moisture 
(from [5]). 
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a) Sandy Loam c) Silty Clay 
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Figure 2.14 Backscattering coefficient at 4.6 GHz, 10 degrees, HH polarization 
as a function of 0-1 cm gravimetric soil moisture (a) sandy loam, 
(b) silty clay loam, and (c) silty clay. Saturated soil conditions 
in the field of silty clay loam are not included in linear regres- 
sion (from [5]). 


expression relating the field capacity to the relative composition 
of the soil (sand and clay contents). 

The impact of soil texture on the sensitivity of a® to moisture 
content was recently investigated by Dobson and Ulaby [5]. Figures 
2.15 and 2.16 depict the observed dependence of on moisture for 
several different sampling depths, with moisture expressed volumetri- 
cally in Figure 2.15 and as percent of 1/3-bar moisture in Figure 2.16. 
Comparison of the two sets of linear regressions clearly demonstrates 
the utility of the conversion to percent of 1 /3-bar moisture for re- 
moving the dependence on soil texture. The overall response to soil 
moisture is summarized in Figure 2.17, where plots of the linear cor- 
relation coefficient between a°(dB) and moisture in the top 5 cm soil 
layer are shown for different soil moisture indicators including grav- 
imetric, volumetric, percent of 1/3-bar moisture, percent of 1-bar 
moisture and percent of field capacity, where field capacity is esti- 
mated on the basis of the expression given by Schmugge [25]. Based on 
the results of the above study, soil moisture given in future sections 
is expressed as percent of field capacity, primarily because the major- 
ity of data acquired in past years does not include soil-tension- 
versus-moisture curves, and therefore Schmugge 's equation is used as 
an approximate estimator of 1 /3-bar moisture. 

2.4 V egetation Cover 

Experimental observations indicate that for angles between nadir 
and about 20®, vegetation cover exercises a minor influence on o° for 
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Backscatter Coefficient (dB) Backscatter Coefficient o® WB) 


ORIGINAL HU"' r-j 
OF POOR QUALITY* 


Figure 

V . ’ ■ 



Sandy Loam 

Silty Clay Loam 

Silty Clay 

All Soils Combined 




.16 Linea)i* regression fits of radar response at 4.625 GHz, 10 degrees, 
and HH polarization to percent of 0.33 bar water content in the 
0-1, 0-2. 0-5, and 0-9 cm soil layers. Regression results are 
based on all data obtained for each soil texture and an equally 
distributed combined texture data base (from [5]). 
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Correlation Coefficient (o®. M 


OF POOR QUALITY 


■ Gravimetric, Mg 
o Volumetric, Mv 

A Normalized Moisture, Mnfor Mtat 1/3 bar 
★% of Field Capacity, Mfc 
□ Normalized Moisture, M„ for Mt at 1 bar 


Frequency (GHz) 


Figure 2.17 


Linear correlation coefficient as a function of frequency for 
various 0-5 on soil moisture indicators with a° at 10 degrees 
incidence angle and HH polarization. Soil moisture is expressed 


as gravimetric Mg, volumetric M^, percent of 0.33 bar moisture, 
percent of estimated field capacity M^^* and percent of 1.0 bar 
moisture. Data is from the combined multi texture data base 


(from [5]), 




frequencies below 8 GHz. Figure 2*18 shows a° versus Mp^ for fields 
of corn, wheat, milo and soybeans, observed by the University of Kansas 
MAS system at 4.25 GHz, and Figure 2.19 is a scattergram containing 
324 measurements acquired by the MAS system ov^r a three-year period, 
of which 181 points are for bare fields and 143 points are for vege- 
tation-covered fields. It 1 $ observed that the slope of the bare 
regression line is slightly higher than that of the vegetation regres- 
sion line, which is attributed to the slight attenuation by the 
vegetation cover. 

Further verification of the above "apparent" independence of 
vegetation cover was obtained from airborne scatterometer data. 

Figure 2.20 shows g° versus Mp^ for 29 vegetation-covered fields 
and 21 bare fields as observed by the 4.75 GHz C-130 scatterometer. 

The data include only fields with no row-periodicity and fields with 
row periodicity that were observed in a direction parallel to the 
row direction. 

2 .5 Composite Effects of Scene Parameters 

The 50 fields of Figure 2.20 are treated as a single category 
in Figure 2.21. In addition, separate regression lines are shown for 
wheat stubble fields and for non-wheat fields observed with the look- 
direction orthogonal to the row direction. The reason for separating 
wheat stubble from the other fields is that the row spacing was 30 cm 
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OHIGIKAI PAG2 
OF POOR QUALITY 


Field Capacity as a Function of o® 



Figure 2.18 Soil moisture content in percent of field capacity as a function 
of ao(dB) at 4.25 GHz, 6 = 10® and HH polarization (from [14]). 


ALL: N = 324, r = . 883 , <r"(dB) = . l33Mpj.(0-5) - 14.34 
BARE: N * 181, r = .849, «*“(dB) = .l48M__(0-5) - 15.96 
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Figure 2.20 Radar backscatter response to soil .noisture based on airborne 

measurements at 4,75 GHz, e = 10® and HH polarization (from Cl8]). 
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Figure 2.21 Variation of with percent field capacity for different covers 
at 4.75 GHz, HH polarization and e = 10®, (from [18]). 


and the row depth was 5 - 8 cm for the wheat stubble fields, in contrast 
to the non-wheat fields whose periodicity and row depth were typically 
around 80 cm and 7-12 cm, respectively. As expected, the different 
row structures resulted in different regressions. 

Evaluation of the row-direction effect on the relationship between 
ff° and Mpjj is shown in Figure 2. 22, where angular plots are given of 
three parameters; (a) correlation coefficient between and Mp^, 

(b) sensitivity of to Mpj. (slope of regression line) and (c) Y- 
intercept (value of c°(dB) at Mp^ « 0). Ideally, one would like to 
see the three sets of curves --representing parallel (//) look direction, 
J_ look direction for wheat, and J_ look direction for non-wheat— identical 
for all three parameters. Unfortunately, this is not the ease for HH 
polarization. However, this ideal situation is almost a reality for 
HV polarization (Figure 2.23). 

2,6 Summary of Sensor-Configuration Selection 

Based on the previous sections, the following conclusions are 
reached: 

a) To minimize the effects of vegetation cover on the radar 
response to soil moisture, angles below 30® and preferably below 20® 
should be used. 

b) In the absence of periodic row patterns, the following con- 
figuration appears to provide optimum performance with regard to 
sensing soil moisture content: (a) frequency: 4-5 GHz, (b) angle of 
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Y’Intercept (dB) Sensitivity (dB/%FC) Correlation Coefficient 



Figure 2.22 Correlation coefficient, sensitivity and Y-intercept for linear 
regression of a°(dB) against percent field capacity, at 4.75 GHz, 
HH polarization (from [18]). 
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Figure 2.23 Correlation coefficient, sensitivity and Y-intercept for linear 
regression of a°(dB) against percent field capacity, at 4.75 GHz 
HV polarization (from [18]). 





incidence range: 10-20®, and (c) polarization: HH. On the average, 
the ntoisture being sensed is that contained in the top 5 cm layer, and 
the moisture content is expressed as percent of field capacity of 
that layer. 

c) If periodic patterns are present and if the row-depth-to-period 
ratio is smaller than about 1/10, o° is approximately insensitive to 
look direction for frequencies higher than 4 GHz, and therefore the 
configuration given in (b) may be used, although son« improvement in 
correlation to moisture may be obtained by operating with HV polariza- 
tion instead. The above row-pattern characteristics are typical of 
conditions that prevail during most of the growing season In dry-land 
farming regions. 

d) In regions where irrigation practices are common, row-depth- 
to-period ratios may be as high as 1/4, in which case the variation 
due to look direction becomes very large, thereby introducing ambi- 
guities in the estimated value of soil moisture. In this case, HV 
polarization should be used because of its weaker sensitivity to look 
direction (in comparison to HH polarization). 

3.0 SYSTEM CONFIGURATION STUDY 

3.1 Introduction 

Imaging radars for soil moisture determination may take on several 
different configurations depending upon the resolution needed. A major 
purpose of this study Is to determine the resolution required and, con- 
sequently, its impact on radar system parameter selection. 

Should fine resolutions be required, no alternative to the nearly 
fully focused synthetic aperture radar (SAR) exists. On the other hand, 
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as the required resolution becomes coarser, the options available to the 
system designer Increase. Various forms of partially focused and un- 
focused SAR may be used. When the resolution requirement Is coarse 
enough, one may consider using a real-aperture sidelooking radar (RAR). 

A "RAR" Is much easier to build than Is a "SAR," so this option must be 
examined seriously. Furthermore, when the resolution requirement is 
relatively modest, one may consider a combined microwave radiometer and 
scanning synthetic aperture radar (RADISAR). The resolution for the 
radar can be of the order of hundreds of meters, whereas that for the 
radiometer will be a few kilometers or tens of kilometers, but the SAR 
picture is embedded within the radiometer cell so that the two may be 
used jointly for the soil -moisture determination. 

As the resolution requirement is tigntened, the required power, 
data rate and processing complexity for the radar all increase. The 
power required for the RAR Is usually quite small at the swath-width 
that one uses for soil -moisture determination. The power required for 
a partially focused or unfocused SAR may also be. quite small under these 
conditions. The power for the RADISAR may be extremely small because 
the total antenna area may be larger for a RADISAR than for a sidelook- 
ing radar that does not scan its beam. 

The complexity of the processing for fine-resolution SAR precludes 
on-board processing. However, for coarser resolution SAR, the complex- 
ity Is reduced by such a large amount that one may think very seriously 
about on-board processing, which in turn reduces the required telemetry 
rate to an almost negligible level. Therefore, one must examine the 
resolution requirement extremely carefully and determine from this 
whether the advantages for other purposes of fine resolution should be 
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traded against the advantages of coarse resolution in terms of low 
power, on-board processing, and low telemetry rate. 

In this study a particular spacecraft configuration has been con- 
sidered as shown below: 


Angles of incidence e: 
Spacecraft height h: 

Frequency f: 

Wavelength M 
Noise figure F: 

Signal-to-noise ratio S^(min): 


70 - 22® 
600 km 
4.75 GHz 
6.315 cm 
4 (6 dB) 
4 (6 dB) 
2 (3 dB) 


Loss factor a: 

Scattering coefficient o°(min): 0.008 (-21 dB), @ e - 22® 

Antenna length D: 8.7 m (8.7 m and 15 m for RAR 

and 5 m for RADISAR) 


Antenna radiation efficiency n: 0.75 
These values are considered to be reasonable ones for implementation. 
The value of the scattering coefficient is a bit higher than usually 
applied in radar design, but with the small angles of incidence and 
typical land backscatter observed in our ground-based measurement pro- 
grams, we feel sure that this is an adequate design level. 

From the fundamental geometric parameters one can calculate cer- 
tain other derived geometric parameters that are used repeatedly in 
the discussion below. These include (Figure 3.1); 

Ground swath width S„: 143 km 

9 

Slant swath width R.: 38.5 km 

Antenna pointing angle range <S: 6.39® - 20.02® 

These parameters are common to all of the systems considered 
here, the primary variables being ground resolution in both across- 
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track direction r. and alongtrack direction r.. The number of Inde- 

y « 

pendent samples averaged (number of Independent looks) N Is a very 
important factor that must be considered In evaluating the performance 
of a system of this kind. Section 3.2 discusses this matter. 

3.2 Equivalent Resolutions 

An experiment performed at the University of Kansas some years 
ago [26] indicated that one can define a resolution volume from which 
the Interpretabllity of an image can be determined. The experiment 
involved the use of human Interpreters evaluating Images for features 
of significant geographic importance; therefore, one cannot be assured 
that the same criteria would apply for a measurement system for soil 
moisture. However, it appears likely that something similar should 
be applicable; so, lacking this improved description of the usefulness 
of an image, one can take advantage of the previous study to rate the 
various simulated Images produced in this study. 

The most significant conclusions of the referenced study are: 

(a) Square and rectangular pixels having the same area are 
equally interpretable. 

(b) A gray-scale resolution can be defined such that it can 
be used with the alongtrack and acrosstrack resolutions 
to produce the resolution volume. 

The resolution volume is defined as 

V = = fe (3-1) 

where V is the resolution volume, r^ is the alongtrack resolution, r 

a y 

is the acrosstrack resolution, r is the newly defined gray-scale reso- 
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lution, and is the pixel dimension of a square pixel of photographic 
quality having equivalent interpretability. The referenced study found 
that a numerical measure of the interpretability of an image is related 
to the resolution volume by an exponential; presumably a similar rela- 
tion may apply for the soil moisture problem, but, in fact, some of the 
simulations reported later indicate that the larger resolution volumes 
may be better for the soil moisture problem. 

The gray-scale resolution rg was defined after considerable ex- 
perimentation with other measures as 

p B signal exceeded 10% of time, v^(2N) distribution ( 3 . 2 ) 
g signal exceeded 90^ of time, x^(2N) distribution 


That is, the gray-scale resolution is the ratio of the signal level 
exceeded 10^ of the time (brighter pixels) to the signal level for 
which 10% of the pixels are blacker. Thus it is a ratio of the bright- 
est-to-darkest picture elements in the speckle pattern of the image - 
a sort of dynamic range for the speckle of the image. One may also 
think of it, for the soil -moisture problem, as a measure of measurement 
uncertainty. 

When the number of independent samples exceeds about four, a simple 
expression can be used to approximate r_: 
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r ~ ^ 

g " /IT - 1.282 


(3.3) 


The value of rg allows one to determine more readily the effect of hav- 
ing a given number of independent samples than one can determine either 
from the standard deviation associated with that number of samples 
(standard deviation is not a good measure for a chi-square distribution) 
or from simply looking at the number N. 
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The fact that one my obtain equivalent results with a rectangular 
and a square pixel having the same area means that one may perform trades 
in the design between resolution in the alongtrack and in the across- 
track direction, as long as the area of the pixel remains the same. In 
the referenced study, this was found to hold in farmland and other areas 
up to a ratio of 10:1 between the length and width of the pixel, 

3.3 System Configuration Studies 

The three types of systems considered here are the real -aperture 
radar (RAR), the synthetic-aperture radar (SAR), and the scanning radar- 
radiometer (RADISAR). Since the simplest system is the RAR, it is con- 
sidered first, followed by the SAR, for which the most examples are pre- 
sented, and finally by the RADISAR, for which only one example is given, 

3.3.1 Real -Aperture Radar 

Figure 3.2 shows the configuration of a real -aperture radar for 
use in space. The typical real -aperture radar used in aircraft uses 
a noncoherent pulse transmitter, but the peak power required for such 
a transmitter is usually too high in space applications thus, a co- 
herent with a pulse-compression scheme ("chirp") must be used. A fre- 
quency synthesizer is therefore required to generate coherent local - 
oscillator and transmitter signals. The coherent transmitter signal 
is fed to a chirp generator that produces an expanded-length reduced- 
amplitude pulse. Either frequency-modulation or binary-phase-code 
modulation may 'le used for this purpose. 
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Antenna 




Figure 3.<? 
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This pulse Is then amplified to the required level and fed to the 
antenna through the duplexer that allows use of a single antenna for 
transmitting and receiving. The antenna is long and narrow^ in this 
case 8.7 m or 15 m long by 40 cm high (the height determines the swath 
width of 7® - 22®). The length determines the alongtrack resolution 
and consequently should be as great as is feasible. 

The signal returns to the antenna and is fed through the duplexer 
to a low-noise amplifier, through a mixer and intermediate-frequency 
amplifier, to a de-chirping system that converts the signals from the 
stretched pulse into signals with the appropriate range resolution. 

This signal is then detected and digitized into a sequence of binary 
bytes, each of which corresponds to a given range. 

The sequence of bytes is fed into a shift register and, as each 
successive pulse is received, the bytes corresponding to the same range 
from this and preceding pulses are added by recirculating the signal in 
the register and adding on a range-element by range-element basis at 
the input. This register may collect all the signals from the pulses 
corresponding to a given alongtrack resolution cell with the same 
amplitude weighting and produce one output for that alongtrack cell. 
Ancther configuration recirculates each pulse with a slightly less- 
than-unity gain, thereby resulting in a continuous output, since the 
effect of the earlier pulses gradually dies out. The output is used 
both to feed the telemetry-or recording-system and to feed a gain- 
control computer. 

The gain-control computer may have some reference variation in 
gain with range that is expected to be typiciil of the terrain imaged. 
However, it also uses the output of the recirculating register to pro- 
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duce an average signal return versus distance curve which can then 
be used to control automatically the gain of the amplifier so that 
the average return at each range is the same. Of course, the gain 
actually used for a particular range must also be telemetered to the 
ground so that the actual amplitude involved can be established for 
the measurement of soil moisture. 

The form of the radar equation used to determine the average trans- 
mitter power is 


.■s 

t 


8.kTj,FSu6^h^a 

a°cos^9 


(3.4) 


where: 

k * Boltzmann's constant 
Tq ■ the reference temperature 290 K 
F * the noise figure of the receiver 
S » the required signal-to-noise ratio 
u = the velocity of the spacecraft 

* the vertical beamwidth of the antenna 
h « the height of the spacecraft 
o = the loss allowance 
n * the radiation efficiency of the antenna 
D « the length of the antenna 
X = the wavelength 
■ ry « the acrosstrack resolution 
<7® = the scattering coefficient 
9 = the angle of incidence (relative to vertical) 

In deriving this equation the following assumptions have 
been made: 
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Antenna gain 6 Is 


where is the alongtrack beamwidth. 

Q 

Pulse repetition frequency PRF Is 
PRF » 2u/D 


(3.5a) 


(3.5b) 


6^ » X/D (3.5c) 

■ slant 4so l ut<bn (3.5d) 


Equation (3.4) has been used In determining the required average power 
In the subsequent tabulations. The peak power required depends upon 
the duty cycle which, in turn, depends jointly upon the PRF and the 
amount of pulse compression used. These are Independent choices that 
may be made by the designer. 

The number of Independent samples potentially available for RAR 
and SAR and actually available for the RAR has been shown [27] to be 


N 


a 



(3.6) 


An Important factor for fine resolution Is the rate at which the 
A/D converter must operate to digitize the range resolution cells. This 
Is 


A/D rate * 



c_ 

2rySlne 


min 


where 


(3.7) 
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Tg ■ effective pulse-length In seconds (after compression) 
c • speed of light 

%in " 0^ incidence at the inner edge of the 

swath, in this case V, 

Another important quantity for a spacecraft radar is the required tele- 
metry rate (or recording rate if an on-board recorder is used). This 
Is simply the number of pixels per second and is given by 


telemetry rate 


'•R''a 


(3.8) 


where is the slant swath width, and r^ is the slant range resolution 
(r^ * rySine). 


3.3.2 Synthetic Aperture Radar 

Figure 3.3 illustrates the synthetic aperture radar configuration 
assumed. The radar itself is a common form; in the processing many 
options are available and only one is indicated, since only one needs 
to be shown to illustrate the relative complexity of the different op- 
tions available. 

The signal transmitted is essentially the same as far the RAR shown 
in Figure 3.2, and the receiver is the same through the de-uni rp system, 
with one exception. Slight variations in velocity or pointing angle 
of the antenna can cause Doppler frequency shifts that are a probl^im 
with SAR. These can be compensated for by adjusting the local oscilla- 
tor either at the first mixer or at the product detectors shown for the 
I rjnd Q channels. Hence, a signal from accurate navigation sensors must 
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Figure 3.3 SAR configuration assumed. 
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be fed into the frequency synthesizer. Following the de-chirp system, 
the synthetic-aperture radar differs from the real -aperture radar. The 
signal is beaten down to zero by feeding a local oscillator signal at the 
intermediate frequency into the product detectors I and Q. These repre- 
sent the in-phase and quadrature components necessary to distinguish 
positive Doppler frequencies from negative Doppler frequencies in the 
output. The processor has independent channels from this point on until 
everything is recombined in the unit at the end, called "Temporary Store 
and Multiplexer." At this point the in-phase and quadrature components 
are combined to make a single spectral component corresponding to a 
particular Doppler frequency which, in turn, corresponds to a particular 
angle off the side looking direction. 

Since the outputs from the receiver are in a range sequence for 
each pulse, but SAR processing requires an alongtrack sequence at each 
range, a corner-turning memory is required. This memory is basically 
a matrix in which the signals are entered in range sequence as each 
pulse returns, and the signals for a particular range are read out in 
azimuth sequence. For a focused or partially focused synthetic aper- 
ture radar, these azimuth outputs are then mixed with a slowly varying 
chirp signal that represents the Doppler frequency shift due to motion 
of the target from a positive Doppler frequency when the 
target is in the beam through zero, to a negative Doppler frequency as 
the target leaves the beam. Thus, these mixers in essence amount to 
correlators between the referenced signals and the individual target 
"chirps." The output of each mixer is then Fourier- transformed 
through an FFT computer to produce the equivalent of the output of a 
bank of filters. Each frequency output at this point corresponds to a 
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different alongtrack element contained in the beam. A separate mixer 
and FFT is required for each range element so that the number of FFT 
computers needed may be very large. However, these computers operate 
relatively slowly, so configurations may also be established whereby 
faster FFT computers can be shared among different channels. For 

simplicity in discussing the complexity of the processor, we assume 
that each range channel has its own FFT unit. For an unfocused SAR, 
no reference function is required. 

The outputs of the FFTs are fed to the Temporary Store and Multi- 
plexer that combines the I and Q outputs, as mentioned earlier, 
and also serves to sort out the returns at the different azimuth angles. 

This permits superposition of the different independent looks at a particular 
point on the ground, since these looks are obtained at different times. 

For example, when a target enters the antenna pattern, it is at a max- 
imum forward angle and consequently a maximum positive Doppler frequency. 

In a multilook processor, the signal is only observed for a given look 
during part of the passage of the antenna beam by the target. At the 
end of that time, the look is completed and the signal is placed in 
the store. The next time this target is observed, it is at a different 
azimuth angle and comes out of a different "filter output" of the FFT. 

Hence, the Temporary Store and Multiplexer must account for this, and 
first take the output for a given target from the highest-frequency 
"equivalent filter" of the FFT, then from the next-highest-frequency 
"equivalent filter," and so on until all of the independent looks have 
been properly averaged. Obviously, this is a relatively complicated 
process, so the Temporary Store and Multiplexer is a fairly complex 
computer. 
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It should be emphasized that many other configurations of synthe- 
tic-aperture processor could be considered. However, the one shown 
here is adequate for illustrating the relative complexity of such sys- 
tems. If one simplifies the system by some alternate design, the simp- 
lification will be in proportion for the different alternative radar con- 
ftguriiltions considered here. 

The values calculated for the real-aperature radar for average 
power, number of independent samples available, analog-to-digital pro- 
cessing rate, and telemetry rate use the same relations as those for 
the synthetic aperture, so that Equations (3.4) through (3.8) apply 
to ’the SAR as well as to the RAR. 

The length of the FFT required is simply the number of pulses con- 
tained within a synthetic aperture at the maximum angle of incidence, 
where the synthetic aperture is the longest. This is the ratio of the 
length of the alongtrack resolution cell for the real aperture (along- 
track illuminated area) to the spacing between potential independent 
samples D/2, divided by the number of independent looks into which the 
potential synthetic aperture is subdivided. The resulting expression is 

FFT Length = — (3.9) 

ND^cose 

max 

The number of FFTs involved is simply the number of range cells and 
this is the ratio of the slant swath width to the slant range reso- 
lution rj^. That is 

R;, 

Number of FFT's = — (3.10) 

*^R 
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The size of the corner-turning tnemory, in principle, would be 
simply the product of these two, but because of the time sequence in 
which signals enter the radar, it must be somewhat longer than that. 

A conservative number is twice the product of the quantities described 
in Equations (3.9) and (3,10). Therefore, the size of the corner- 
turning memory is 


Size of C.-T. Memory = 


4xhR. 


2 

0 r„(0^,.„)sine„-._cos0„^ 
y min^ mm max 


(3.11) 


It should be recalled, however, that this corner-turning memory 

and the FFTs must be duplicated for the in-phase and quadrature channels. 
Furthermore, when excess bandwidth is used to obtain additional inde- 
pendent samples beyond those available in the alongtrack direction, 
the number of corner-turn memories and FFTs is multiplied by the excess- 
bandwidth ratio needed to achieve the required number of independent 
samples. Therefore, the processor in such a situation becomes extremely 
large and complex. 

The total number of independent samples is the product of the 
number obtained by excess range bandwidth and the number ob- 
tained in the alongtrack direction 

N = 

The number required in the range direction is the excess-bandwidth 
ratio mentioned above, and N is given by Equation (3.6). 

di 

In Figure 3.3, reference functions are shown to be required for 
the SAR processor. The number of different reference functions re- 
quired depends upon the depth of focus, which is given by 
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(3.12) 


Depth of Focus ' 2r;/X 

a 

The depth of focus is defined in the slant-range direction so that the 
number of reference functions required is simply the ratio of the slant 
swath width to the depth of focus or 

\R_ 

Number of Ref. Funcs. * — ^ (3.13) 

For the first three cases discussed in the next section, more than one 
reference function is required. For the remaining cases of partially 
focused synthetic aperture a single reference function suffices, and 
for the unfocused synthetic-aperture, no reference function is required. 

3.3.3 Scanning Radiometer-SAR 

Figure 3.4 illustrates the scanning synthetic-aperture configura- 
tion that can be combined with a radiometer to produce the RADISAR. 
Basically, this system uses a scanned-array antenna and filters to 
separate the SAR signals from the radiometer signals; beyond that it is 
the same (except in some details of processing) as the other SARs and 
any radiometer design may be chosen. A larger antenna area may be used 
for the RADISAR than for the other configurations discussed here be- 
cause of the fact that the antenna scans from one position to another, 
so the antenna height is not constrained by the swath. For this reason 
the power required is less. 

3.4 System Tradeoffs 

Various systems have been postulated for use in the soil moisture 
problem, with the reference system being a fully-focused synthetic- 
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Figure 3.4 RADISAR configuration. 


57 



aperture radar with antenna 8.7 meters long and with a square pixel 
at 7®, Table 3.1 summarizes some of the characteristics of the Systems 
considered, and Table 3,2 gives more details. The characteristics 
shown in Table 3.1 are: the equivalent square photographic quality 
pixel dimension r^ discussed in Section 3.2, the diirensiorts of the 
radar pixel at 22®, the number of looks used in the postulated system, 
the required average power, the speed required for each analot,-to- 
digital converter, and the telemetry rate required if all processing 
is done on-board. If processing is not done on-board the telemetry 
rate is much higher (except for Case 1). 

Case 1 is the fully focused reference case. Mote that to achieve 
the 4.4 m resolution at 7° in the range direction, a 1.4 m resolution 
is needed at 22®. The power required is quite large and the A/0 
converter rate and telemetry rate are exceedingly large. 

In Case 2, a partially focused SAR with 10 x 10 m resolution at 
22® is compared with Case 1. Since it uses 12 looks to reduce the 
variance of the signal, and has to obtain many of those looks by ex- 
cess bandwidth, its power requirement is also large, but the A/D 
converters operate at a rate about one order of magnitude less than for 
Case 1, and the telemetry rate is also an order of magnitude less. The 
equivalentphotographic-quality pixel dimension is not much more for 
Case 2 than for Case 1. This kind of resolution is probably entirely 
too fine for a soil -moisture radar. 

Case 3 shows the resolution likely to be the finest that might 
be considered for a soil-moisture radar, and it has a 30 x 30 m resolu- 
tion at 22°. Because it requires excess bandwidth, it also has a fairly 
high power (800 W), but its telemetry rate is down significantly. 


58 



In Case 4, the processing is greatly simplified, but the range 
resolution is kept the same. Because the excess bandwidth is no^ re- 
quired for Case 4, and is for Case 3, the power is less in Case 

Case 5 is a partially-focused SAR with 100 x 100 m resolution r ’ 

22® and 23 independent samples. This might be a realistic design. How- 
ever, as indicated by the resolution study, this kind of resolution is 
probably not required; hence, the unfocused cases. Cases 6-8 and 11, are also 
shown. Note that the po',er required can be very small for these cases. 

Case 8a is shown in Table 3.1 to illustrate the effect on the power of 
exchanging the azimuth resolution for the range resolution. Of course, 
this is done at 22° and at 7® it is not quite as good a system. How- 
ever, the power is reduced from 46 to 13.8 watts. 

Four RAR cases are shown, two with 8,7 m antenna lengths like the 
SAR, and two with 15 m antenna lengths. Cases 9 and 10 are designed 
to achieve a 7° equivalent photographic pixel dimension of about 1 km. 

As a result, the pixels have a very high length-to-width ratio, prob- 
ably higher (at least in Case 9) than one can use to extend the results 
of the study indicated in Section 3.2. The 20:1 ratio in Case 10, how- 
ever, may be reasonable. The other two cases have arbitrarily self^cted 
rather large range resolutions at 22°, so that the effective square- 
photographic-quality-pixel dimensions for these cases are very large. 

Note, however, that the powers are quite small. 

Case 14 7s the only case of the RADISAR considered. For the RADISAR, 
the 40 cm height of the other antennas needed to achieve the 7° - 22° 
swath width has been replaced by a 5 m height, and the antenna has been 
shortened to 5 m. This square array will produce nearly circular radio- 
meter patterns on the ground. With this size antenna, 19 scan positions 
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are required for the scanning synthetic-aperture radar, which means that 
there are 19 radiometer measurements made across the swath, Because of 
the larger antenna, the average power is only 0.66 watts. 

Table 3.2 gives more details of the systems described in Table 3.1, 
and also repeats some of the significant factor« listed in Table 3.1. 

Table 3.3 gives additional information regarding the complexity 
of the synthetic-aperture processors, namely the length of and number 
of FFTs and the size of the corner-turn memory (in bytes). Note that, 
in each case, a number for the FFT length that is a multiple of two 
is given. This is the actual length of the FFT, although some other 
transforms do not require that the length of the input be a multiple 
of two. In Case 1, both 1024 and 2048 are listed because the 1080 is 
so close to 1024 that only a small compromise would be made in using 
it. The corner^turn memory need not contain the additional zeros that 
enter the FFT. 

A column is also i ^eluded in Table 3.3 listing the number of FFT 
systems and corner-turn memories required. Thus, the numbers in columns 
3 through 5 are for each set of these items and the total number re- 
quired is that indicated in the last column. 

The tradeoffs between resolution and the equipment required for 
the system, as well as the telemetry rate, are best illustrated by 
comparing the values of equivalent pictorial quality resolution with 
the other quantities. Figure 3.5 illustrates the relation between this 
equivalent pixel dimension and the power and telemetry requirements. 

For the SAR, the power requirement decreases rapidly as the pixel dimen- 
sion increases. Since the power required for the SAR depends primarily 
on the range resolution, cases involving similar pixel areas such as 
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8 and 8a have similar photographic pixel dimensions but different 
power requirements. On the other hand, cases where range resolution 
applies to two different alongtrack resolutions, such as 3 and 4jOr 5 
and 6, require the same power (except for the increase in 3 due to use 
of excess bandwidth). 

Although RARs follow the trends of the SAR reasonably well, the 
power requirements can differ quite significantly since the RAR power 
requirement also depends upon the range resolution. For instance, to 
achieve an r„ of about 1000, Case 9 requires a 71 m across track resol u- 
tion at 22®, which causes the average power to go up to 194 W, a value 
much higher than required for the same equivalent pixel dimension with 
the SARs of Cases 8 and 8a. The RADISAR is well off the trend because 
its antenna area is not constrained to a reasonable length multiplied 
by a fixed height. Clearly the RADISAR is the least power-consuming of 
the options, but the processing for comparable r_ is more complex in ways 
not indicated by this kind of comparison. It also requires a larger 
and more complex antenna. However, this antenna makes possible the 
use of both lower power for the radar and of the radiometer along with 
the radar, which may be a significant advantage. 

3.5 Relation of Tradeoffs to Simulations 

Not all of the examples illustrated in Section 3.4 could be simu- 
lated. The simulation process is an expensive and time-consuming one, 
and the simulations must start off with a spatial resolution for the 
data base of 20m x 20m. For this reason, an actual simulation of Case 1 
was impossible. An actual simulation at 7® of Case 2 was not possible, 
but one with a comparable equivalent square-photographic-quality pixel 
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was possible. However, the same could not be said about Case 2 at 22®, 
since the effective resolution there was finer than could be readily 
achieved with the data base. Also, although the techniques are available 
for degrading the image from a 20 x 20 image to any other value with 
larger pixel dimensions and with any number of equivalent independent 
samples, some methods are easier to carry out than others, and this 
tends to constrain the number of independent samples in the simulations. 

Two sets of simulations were performed, covering the angular ranges 
of 7.5®-9.3° and 11°-12.8®. A third set, covering the 17.96®-19.6° range 
was planned but was not performed due to time and cost constraints. In 
Table 3.4, a comparison is made between the simulations performed and the 
nearest equivalent system-study cases. 

The simulations labeled A all are at the spatial resolution cor- 
responding to the data base. They are all approximately equivalent in 
r to Case 2 at 7®. That is, the value of r„ obtained with the simula- 

c e 

tion was not significantly poorer (within 20^) than the for Case 2 
at 7®. Note that .simulation B is approximately equivalent to Case 4 and 
simulation C is approximately equivalent to Case 8. 

The soil -moisture-measurement implications of the simulations 
are discussed in other sections. However, one c^ say that simulation 
B has a resolution adequate to distinguish fairly easily the features 
of the image, whereas simulation C makes many of the features indis- 
tinguishable, although it is still quite useful for soil-moisture meas- 
urement. Simulation B has a considerably better resolution than the 
RADISAR case studied or any of the RAR cases studied. Simulation Case C, 
although most directly comparable with Cases 8 and 8a, is also reason- 
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TABLE 3.4 

Simulations Compared with Nearest- System Case 


Simulation 

Case 
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Angle 
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(m) 

N 

>^e 
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fg at 6 
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CM 
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GO 

20 

12 

29.5 

2 

25 at 7® 

A2 

n.l®-12.8® 

12.1® 

20 

12 

29.5 

2 

25 at 7® 

B1 

7.5®- 9.3® 

8.2® 

90 

23 

125 

4 

122 at 7® 

B2 

n.l®-12.8® 

12.1® 

93 

27 

96 

4 

94 at 12® 

Cl 

7.5®- 9.3® 

O 

CM 

• 

00 

900 

2250 

974 

8 

976 at 7.5® 

C2 

n.l®-12.8® 

12.1® 

900 

2800 

753 

8 

750 at 13® 
















































ably like Cases 9 and 10 for the RAR and gives a larger equivalent 
pixel dimension than Case 14 (the RADISAR). 

Thus, one can state that the simulations conducted are reasonably 
representative of a fine-resolution case (2), a medium-resolution case 
(4), and a relatively poor-resolution case (8). Insufficient time 
and resources were available to simulate the still -coarser resolutions 
associated with Cases 11, 12, and 13; these may be quite useful for 
soil -moisture determination, but rather poor for other imaging purposes. 
This kind of simulation should be performed for these poorer resolutions 
and, indeed, extended to the larger cells that are likely to be found 
with a radiometer. Similarly, simulations to show the effect of pixels 
that are rectangular rather than square should be conducted for ratios 
of r^ to ry considerably larger than those used in the simulations con- 
ducted to date. Not only may a maximum value of this ratio exist be- 
yond which the concept of the gray-scale resolution and the resolution 
volume does not apply, but some maximum value for the alongtrack or 
acrosstrack pixel dimension may exist beyond which the results will be 
of little value. 

4.0 SIMULATION STUDY 

4.1 Radar Image Simulation 

Image simulation has been used as a research tool to help deter- 
mine, sensor parameters for an orbital altitude radar whose images are 
to provide soil moisture information. The sensor parameters include, 
but are not limited to, operating frequency, polarization, angle of 
incidence, resolution, averaging, and data format. The investigation 
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of these parameters and the ability of spaceborne radar to provide 
quantitative soil moisture information was carried out by varying both 
sensor parameters (angle of incidence, resolution, averaging) and ground 
conditions (varying amounts of soil moisture). The frequency and polar- 
ization for the modeled sensor were chosen according to studies of ter- 
rain backscatter reviewed earlier in this report. 

The following sections describe the physical processes of radar 
image formation, and the design and implementation of algorithms 
to generate radar images. Later, the results (images) are analyzed 
quantitatively to determine how well one can estimate soil moisture 
from a processed radar image. A processed radar image is one in which 
known sensor-dependent effects have been removed. For example, antenna 
gain variations in the range dimension are removed by scaling the image 
data according to an angularly dependent function. 

The radar image simulation algorithms used in this study were 
developed by Stiles et al. [28]. fne simulation procedure follows 
closely the actual operation of synthetic aperture radar imaging; thus, 
it is relevant to present an overview of the image formation processes 
herein. Sensor, terrain, propagation channel, geometric, and data pro- 
cessing effects will be treated in Sections 4.1.1 through 4.1.5. The 
overview and simulation will be presented in terms of generating fully 
focused synthetic aperture radar images. Changes necessary to accom- 
modate the production of real aperture and partially focused SAR images 
will be discussed in Section 4.1.5. 
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4.1.1 An Overview of the Radar* Image Formation Process 


The orbital platform for a coherent radar system has known velocity, 
altitude, and attitude; the radar antenna Is usually oriented with Its 
electrical boresight perpendicular to the spacecraft velocity vector. 

The sensor transmits a waveform with known amplitude and phase modula- 
tions which are used In the information decoding process. 

The pulse repetition frequency (PRF), is designed to simultaneous! v 
satisfy range ambiguity limitations (which set the maximum PRF) and 
azimuth ambiguity limitations (which define the minimum PRF) [35]. The az1 
muth dimension is parallel to the platform velocity vector while the 
range dimension is orthogonal to the velocity vector in the most common 
configuration. The power density transmitted into space by means of a 
directional antenna is large in the direction of the antenna main lobe 
and small in the direction of the side lobes. Due to diffraction the 
power density that is incident upon the ground is r'‘'‘iced from the trans- 
mitted value. The scattering and reflection processes which occur next 
direct some portion of the power in the "backscatter" direction, i.e., 
toward the radar receiving aperture. The backscatter cross-section of 
the illuminated terrain (a) incorporates all ground effects needed to 
estimate the average received power by the definition a = 
the ground effects can be further broken into a backscatter coefficient a° 
and an effective ground area term A^^^. Various types of terrain can be 
represented as having different backscatter responses, thus, the action 
of the terrain upon the transmitted pulse is characterized by amplitude 
modulations. Additionally, phase modulations are incurred because of 
the relative motion of the antenna with respect to the ground scatterers. 
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These modulations have been theoretically modeled, and are also used 
in the signal decoding process. 

Each transmitted pulse illuminates a swath of ground; the corres- 
ponding received pulse is mixed with a stored signal to put the informa- 
tion on a carrier of lower frequency. Then the packet of information 
is written as a modulated bipolar video frequency onto a "signal film," 
or it is stored digitally as a matrix row. Each subsequent received 
pulse is written alongside the previous one in synchronism with the space- 
craft's velocity. To relate the fngnal film coordinates to the ground 
geometry, distance across the width of the film represents increasing 
range from the antenna, while distance along the film represents azimuth. 

The processing of the raw data is simple in concept, but complicated 
in practice. It is appropriate to think of the raw data as 
a two-dimensional complex signal having amplitude and phases additionally, 
the signal film consitutes a Fresnel zon'i plate whose diffraction prop- 
erties are used to produce the final image. Two basic operations are 
performed on the data in the orthogonal dimensions of range and azimuth: 
range compression and azimuth focusing. Both are implementable by fre- 
quency-dependent phase shifts; the end-product following these operations 
is a "complex image" referred to as a radar hologram. Either the real 
amplitude or real magnitude of the hologram is usually used to develop 
a film negative for the final radar image.’ 

The resultant radar image has properties much different from those 
of noncoherent imagery. The phase coherence of the sensor not only 
allows fine resolution images to be generated, but it also inherently 
causes a strong correlation to be induced upon the backscattered signals, 
thus making the image representation of the fading signal appear speckled 
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and worn\y. Additionally, the geometry of radar gives it unique prop- 
erties which one must understand for interpretation purposes. These 
geometrical factors will be discussed in the following sections. 

4.1.2 Target Geometry Considerations 

Four geometric effects are of considerable importance in radar 
imaging and the simulation of imagesi foreshortening, layover, shadow, 
and local slope. Additionally, earth curvature can be important, and 
is discussed in Section 4.1.?. Foreshortening refers to object com- 
pression seen in the radar irjftge, and it occurs for positive slopes 
(landform rising along the increasing range direction). Layover is a 
more severe case which is commonly demonstrated by the "flagpole" ex- 
ample in which the incidence angle of the radar and the height of the 
pole dictate that the top of the flagpole is "seen" by the radar before 
the base of the pole. Shadow is the absence of illumination due to the 
presence of trees, buildings, mountains, etc., which occlude the terrain 
that they precede in range. These geometry-related factors are incorp- 
orated in simulation through range calculations with respect to the 
radar, and by neighborhood comparisons. Specifically, average return 
powers computed by the radar equation [29] are apportioned to the cor- 
rect range gates or bins. 

The local angle of incidence, e , is defined as the angle between 
a local normal for a scattering element and the range vector running 
from the antenna to the element. Implicitly, a planar facet exists to 
help defiisS the local normal. The Siimulation algorithms examine the 
terrain and define slopes in the range and azimuth dimensions to deter- 
mine the facet orientation. When 8^ is zero degrees, i.e., the scatter- 
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ing element or cell is oriented perpendicularly with respect to the 
range vector, then the backscatter power is predicted by the known back- 
scatter response for that target class, e.g., pasture, for an angle of 
incidence, e, of zero degrees. The assumption is thus made that the 
backscatter response (a° vs. e) should be examined at to Incorporate 
local slope. Thus, the simulation algorithms use a°(0.) to predict 
return power (through the radar equation [29]). 

Multipath effects can be significant for certain terrain conditions, 
but are not included in the simulation algorithms. Multipath reflec- 
tions from water, followed by scattering from the terrain can cause the 
"far shore" of lakes, rivers, etc., to be brightened. Since our appli- 
cation for simulation is to test soil moisture prediction algorithms, 
primarily for agricultural regions, water-body boundaries are not 
of great interest. 


4.1.3 Earth Curvature Considerations 


The sensor platform altitude, being significant in relation to the 
radius of the Earth, causes the angle of incidence, 0 , measured between 
a normal to the ground within the swath and the range vector, to be 
larger than the radar incidence angle measured between nadir (relative 
to the sensor) and the range vector. As shown in Figure 3.1, e > 6. 

The relationship between these angles of incidence is given by 


0 = sin 



(4.1) 


Because 0 also affects the local angle of incidence, 0 , the Earth's 
curvature is quite significant. The mid-range 5 values employed In 
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the simulations, i.e., 8.4® and 11.9®, correspond to e's of 9.1® and 
13.0®, respectively. 

The geometry effects defined above are modeled and incorporated 
as deterministic influences on the amplitude structure of the reflected 
pulse of microwave energy. They affect the average power calculations 
for eich picture element as they define where the power from scatter- 
ing rsgions belong in the rad:-.r image. The scattering action of the 
terrain, which perturbs the amplitude and phase of the transmitted 
pulse, is, however, modeled as a random phenomenon commonly known as 
fading. For most distributed targets the random point process is 
reasonably described by the Rayleigh distribution [30]. The requirements 
for this case are that no one scattering center dominates the reflec- 
tion/scattering process, that there are many scattering centers within 
the smallest discriminable terrain element, and that the phases of the 
returns from the various centers are Independent [30], This topic will 
be discussed further in Section 4.1.4. 

4.1.4 Simulation of Fading and Other System Effects 

The >"iadar system transmits a narrowband signal and receives a 
scattered signal that is characterized as a narrowband Gaussian random 
process [38], This model for the real signal is consistent with the 
Rayleigh scattering model [31]. That is, the statistical nature of the 
envelope can be generated either by drawing samples from a Rayleigh 
random number generator or by generating and processing two independent 
random Gaussian numbers for the N = 1 case (thus, mimicking the quad- 
rature reception of the SAR signal.) The detected voltage statistics 
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are governed by a negative exponential distribution for N » 1. For 
most of the simulations prepared for this study, the following model 
[32,33] was used to compute instantaneous values of the scattered power 



where Pj^ is the power calculated for each facet on the ground, Fp is 
the average power calculated by the radar equation, N is the number of 
independent samples, and y is a random variable which is distributed 
as an exponential variate or as a chi-square variate with 2N degrees 
of freedom. Thus, if an N * 1 image is desired, two free squares enter 
the calculation, i.e., 

4 ’ Zi + Zg 

where and Zg are independent, normally distributed random numbers, 
p(0,l) (zero mean, unit variance). The larger the number of looks de- 
sired, the higher the probability becomes that the actual received power 
is close to Fp, which is given by [29, 34] as 


_ - Sj(0»<^)Gp(0,(j)) X a°A 


’’r' Pt 




(4.4) 


The assumption was made that the factor 

. Gj(e»i>)Gp(0,i>) 2 


(4.5) 


could be written as a constant, and that furthermore this constant 
value c could be calibrated out of the resultant SAR image by a scaling 
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constant. If a linear system model is applicable for the receiver and 
processor, then this scaling is plausible. Since the terrain surface 
slope Is generally unknown, it is not possible to remove the range R 
and resolution area A terms for estimating 

As discussed earlier, the value of a® used for computing the 
average received power from a resolution cell Is where 

» cos"^(p • n) (4.6) 


p ■ the unit vector pointing from the center of the resolution 
cell to the antenna center 

A 

n * local normal unit vector constructed for a facet model of 
the resolution cell 

Average power was actually calculated for the finest resolution possible, 
approximately 20 meters by 20 meters. That is, for each entry In the 
ground-truth data base, the local angle of Incidence was determined from 
the given elevations of the cell in question and its nearest four neigh- 
bors and from the given sensor position. The corresponding value of 
the backscatter coefficient is drawn from a stored table of o° versus 
0j^. The average received power that Is utilized for each location 1s 


P 


R 



Aa®(e^) 


(4.7) 


where c » 1 and A is defined as 


A 



COSp 


sin(9-'l») 


where 


(4.8) 
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d-- » data base resolution In the azimuth dimension 

aZ 

d^ * data base resolution In the range dimension 
0 ^ angle of Incidence relative to a flat surface 

r 

(not the local angle) 

p ■ angle of the slope of the cell in the plane 
orthogonal to the plane of Incidence 
Ip ■ angle of the slope In the plane of Incidence 

The simulations that were produced are In ground range fonnat, 
and this basically means that range gates were explicitly set up in 
non-equal Intervals; In fact, the range gates were purposely made 
smaller at the near- range edge of the Image. The gate length Is given 
by 



• sine 


(4.9) 


Suppose that the slant range resolution Is 25 meters. Then the range 
gate defined by Equation 4.9 has length 8.6 meters for a nominal angle 
of Incidence of 20°. The purpose of the Inclusion of the sine factor 
Is to compensate for the inherent poorer resolution (in range) that 
occurs In the near range of the scene compared with the far range. 

The local angle of incidence, and subsequently the area and back- 
scatter coefficients, are computed for each terrain cell. The average 
power is distributed among range bins according to the position of the 
terrain facet center relative to the start- and stop-points of the ap- 
propriate range bins. For example, if the range to a certain ground 
facet puts it three-tenths of a range bin away from the start of the 
(k + 1) range bin, then three-tenths of the power for that ground 
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facet win be allotted to the k bln and seven-tenths Is contributed 
to range bln k + 1. 

The power 1n each of the ground-range gates Is used as the basis 
for the calculation of the fading value as In Equation 4.2. As was in- 
dicated in Table 3.4, there were various numbers of looks required; 

In cases where N is greater than 12, the approach via Equation 4.2 was 
used to generate N = 1, 2, or 3 looks. These Images were subsequently 
smoothed to trade resolution for averaging via the application of a low- 
pass filter [28] which implements the scanning mixed integration for SAR 
S/N improvement [36, 37] to achieve the looks calculated for Bl, B2, Cl, and 
C2 (Table 3.4). In the case for which N was required to be 12 looks, 
a Gaussian approximation to the operation of Equation 4.2 was used, 
that Is [39] 




V • ii ' 


(4.10) 


where Z is distributed as a Gaussian random variable with zero mean, 
unit variance, i.e., p(0,l). 

Noncoherent spatial domain smoothing was applied to several of the 
images generated via Equation 4.2. Coherent processing will be defined 
as any operation on amplitude and phase of a radar signal; incoherent 
processing will be defined as operations on amplitudes only or on 
power only. Thus, azimuth focusing and range compression are coherent 
operations; signal processing to improve the signal -to-noise ratio is 
usually done coherently with an Incoherent sum of intensity data at the 
end. 
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4.1.5 Noncoherent Processing to Increase Averaging 


When the number of looks required (see Table 3.4) was large (>12), 

It was Infeas^/ble to generate chi-square random numbers to satisfy our 
needs. However, If we first generate N « 1, 2, or 3 simulations, 20 m x 20 m 
resolution, we would be able to use these as bases for subse- 
quent resolution degradations. Therefore, we applied a spatial domain 
convolution (reported 1n [28]) as suggested by Zelenka [37]. The 
low pass filter derived and applied to N ■ 1, 2, or 3 Images accomplishes 
the scanning mixed Integration operation for speckle reduction. 

Both the scanning and discrete mixed integrators [37] use rion- 
redundant spectral Information for a radar scene and they make Intensity 
Images from the Inverse transforms of several subsections of the spectra. 

A moving window Is applied for scanning, and fixed windows at several 
locations are used In the discrete mixed-integration case. Convolution 
simultaneously degrades resolution while averaging out the variations In 
Image intensity due to coherent fading. Thus, it was possible to In- 
crease the number of looks and to degrade resolution by approximately 
Integral multiples to the values specified in Table 3.4. 

4.1.6 Radar System Parameters Modeled 

Three different radar systems were modeled in this study; these 
consisted of two partially focused SAR's and one unfocused SAR (see 
Table 3.4). Since the ground truth data base pixel spacing is 20 nfieters 
by 20 meters, the finest resolution simulations were generated for this 
resoluti^. The frequency and polarization of the radar were 4.75 GHz; 
horizontally transmitting, horizontally receiving, respectively. 
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Due also to data base size limitations (approxirnately 19,2 kilo- 
meters in range by 16.4 kilometers along the azimuth dimension) it \^as 
not possible to represent the desired swath (143 km along the range 
dimension) by one simulated radar image. It was necessary to use the 
same data base for several different angles of incidence as samples of 
the desired swath (see Table 3.4)'. 

^•2 Determination of Mea n Backscattering Coefficient of Target 
I^o^tibns as of Incidence Angle and lloisture 

The data reviewed in section 2.0 were compiled and evaluated to 
determine algorithms which estimate mean backscattering coefficient at 
approximately 4.75 GHz and HH polarization as a function of target 
class and incidence angle. The general form of the estimation algo- 
rithm is given as; 

- f(e^) + 9 ( 04 ) K Mpg (4.11) 

where 

5® » the mean a® in dB of a given target class 

f(Qj^) * a function of incidence angle 0 which includes row- 
direction effects in certain agricultural crops 

g(6^) « the radar sensitivity to moisture for a given target- 
class as a function of incidence angle 

Mpc ■ the 0-5 cm soil moisture expressed as a percent of 
the 1/3-bar water content 

» the local angle of Incidence 
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Source data at center frequencies between 4.25 and 4,9 GHz with 
HH polarization representing each target class in the simulation were 
used to determine f(o) and g(9) in Equation 4.11. For a given target 
class, values of f(o) and g(o) were established by least-squares linear 
regression of as a function of soil moisture MpQ at each available 

angle of incidence. The linear regression coefficients for bare soil 
and crop classes are given in Table 4.1 at incidence angles of 0®, 10®, 
20® and 30®. A representative scatterplot is shown in Figure 4.1 for 
smooth bare soil at 10® incidence angle. The linear regression co- 
efficients were fitted with third-order polynomials as a function of 
incidence angle to produce f(Q) and 9 ( 0 ) as given in Table 4.2. Assum- 
ing to be continuous and correlated as a function of incidence angle, 
the values in Table 4.2 permit the estimation of mean radar backscat- 
tering coefficient for 0® S 0 s 30® and 0 S Mpj. ^ 150. 

All moisture dependent algorithms incorporate percent of the 
1/3-bar water content Mp^ in the 0-5 cm layer as the pertinent soil 
moisture indicator. The 1 /3-bar water content approximates field 
capacity for certain soils. For the data sources listed in Table 4.1, 
the 1 /3-bar water content of the soil was either measured directly 
[42,443 or estimated from soil textural components by [25]; 

FC - 25.1 - 0.21 (% sand) -f 0.22 (« clay) (4.12) 

where 

FC *• estimated water content at 1/3 bar 
and 

Mpj, * 100 Mg/FC 

where 
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TABU 4.1 

Linear Regression Coefficients, Correlation Coefficients, 
Sample Size, and Data Sources of 4.75 GHz, 

HH Radar Response to 0-5 cm Soil Moisture 
of Target Classes at Selected 
Incidence Angles 





Regression 


Number 

Number 



Sample 

Size 

Incidence 

Coefficients 

Linear 

of 

of Exper- 
imental 





Correlation 

Soil 


Data 

Target Class 

f(e) 

g(e) 


Angle 

Coefficient 

Textures 

Fields 

Sources 

$are Soil, Smooth 

46 

0 

- 5.13 

0.182 

.49 

2 

3 


140,41] 


106 

10 

-17.53 

0.160 

.85 

5 

6 


l40,41, 

. 

106 


-21.01 

0.131 

.83 

5 

6 


40.41. 

Bare Soil, Medium 

46 

49 

wm 

-23.82 

-11.69 

0.112 

0.137 

.66 

.73 

2 

2 

3 

3 



Rough 

49 

10 

-15.49 

0.149 

.89 

2 

3 




49 

20 

-17.47 

•0.128 

.82 

2 

3 



Bare Soil, Rough 

49 

30 

-18.85 

0.114 

.83 

2 

3 



29 

0 

-15.09 

0.157 

.83 

2 

2 


7 itt 


29 

10 

-14.18 

0.139 

.88 

2 

2 




29 

20 

-17.06 

0.145 

.89 

2 

2 


C Ik 

Corn 

29 

30 

-19.82 


.88 

2 

2 



33 

0 

- 7.77 

0.128 

.68 

2 

3 


2 ' 


32 

10 

-11.36 

0.104 

.87 

2 

3 


1 CfF ' 


32 

20 

-14.22 

0.076 

.74 

2 

3 


'43,44] 


33 

30 

-15.57 

0.079 

.69 

2 

3 


‘43,44] 

Milo 

22 


- 9.74 

0.124 

.86 

1 

1 


‘43] 


22 


-12.12 

0.086 

.73 

1 

1 


'431 


22 


-13.48 

0.067 

.54 

1 

1 


43’ 

Soybeans 

22 

BB 

-14.47 

0.062 

.57 

1 

1 


43] 

55 

■EH 

-10.00 

0,181 

.91 

2 

3 


‘43,44] 


55 

10 

-13.60 

0.126 

.89 

2 

3 


‘43,44] 


55 

20 

-14.63 

0.093 

.90 

2 

3 


43,44] 


55 

30 

-16.13 

0.096 

.89 

2 

3 


43,44] 

Wheat 

34 

0 

- 1.68 

0.108 

.72 

2 

5 


'43 ',44] 


34 

10 

-16.02 

0.160 

.91 

2 

5 


‘43,44] 


34 

20 

-12.84 

0.025 

.35 

2 

5 


'43,44] 


34 

30 

-14.20 

0.019 

.23 

2 

5 


[43,44] 
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TABLE 4.2 (continued) 


B. Targets Modeled with no Dependence on Soil Moisture, where 
g(e) « 0 and f(e) i constant. 


Target Class 

f(#)* 

Water Bodies 
Deciduous Trees 

22.82 - 5.126e + 0.2370*- 3,973 X 10‘^e^ 
10 log do’- X cose) 


* f(e) for these targets Is valid for 0° - e - 30®. 


C. Targets Modeled with no Dependence on Soil Moisture or Incidence 
Angle* where g(e) * 0 and f(e) ■ constant. 


Target Class 

Constant Value (dB) 

Railroads 

10.0 

Bridges 

10.0 

Buildings 

10.0 
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Mg « gravimetric soil moisture 


The use of In Equation 4.11 serves to minimize the dependence 
of o® on soil texture at a given gravimetric or volumetric soil moisture. 
In Figure 4,1, each symbol represents a different soil texture for bare 
soil with RMS roughness ~ 2.0 cm. The result of the net linear regres- 
sion result shown for 0 ■ 10®, 


* -17.53 + 0.16 Mp^. (4.14) 

Is not significantly different from the result obtained for any single 
soil texture regressed alone. 

In addition, only data for moisture conditions less than Mp^ * 150 
were used in algorithm development since Mp^ greater than 150 represents 
flooded soil conditions. When moisture conditions exceed saturation, 
the linear algorithm model given by Equation 4.11 for a given incidence 
angle becomes invalid at a given angle since the electromagnetic behavior 
of the target becomes similar to that of a smooth water surface. 

Other target classes show no dependence upon soil moisture in 
establishing a mean response. Deciduous trees, roads, and water bodies 
are modeled with a dependence only on Incidence angle and hence g(e) * 0. 
Cultural targets are modeled with a constant a° for all incidence angles 
and soil moisture conditions. 


4.2.1 Bare Soil Algorithms 

The generalized bare soil algorithms given in Table 4.1 are for 
three classes of surface roughness conditions. For the bare soil case, 
surface roughness is considered for randomly distributed soil aggregates 
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on a planar surface where roughness Is defined by the RMS height In cm 
of the aggregates. For convenience, the range of probable RMS values 
within the simulation test site was subdivided into three roughness 
classes which correspond to distinct phases of soil tillage: 

1) 0 < RMS height - ^,0 cm for smooth bare soil, 

2) 2.0 cm < RMS height - 4.0 cm for medium rough bare 
soil, and 

3) 4.0 cm < RMS height for rough bare soil. 

Within the simulated area, a smooth soil surface would exist for fal- 
lowed fields, especially wheat fields, and would be most prevalent just 
prior to seed bed preparation in the spring. A smooth surface also de- 
scribes that usually existent under certain crop canopies: pasture, 
alfalfa, and wheat. For the simulated area, the rough bare soil condi> 
tion with RMS height exceeding 4.0 rarely occurs under normal circum- 
stances and then only as an immediate consequence of a tillage opera- 
tion such as plowing or cultivation. The rough bare soil condition tends 
to exist only until a precipitation event which reduces RMS height to the 
medium rough class. Thus RMS height > 4.0 cm is transitory and tends to 
exist only in association with dry to moist soil conditions. The medium 
rough soil condition describes the soil surface most prevalent within 
the simulated area. This condition is typical of the soil surface 
beneath the canopies of most row crops: corn, milo, and soybeans. 

A variety of soil textures also are included in the development of 
the bare soil algorithms. The choice of Mp^ as the soil moisture descrip- 
tor effectively minimizes a° estimate error due to soil texture effects. 
Table 4.3 lists RMS values and soil textures of the experimental data 
used to develop the bare soil algorithms given in Tables 4.1 and 4.2. 


TABLE 4.3 

RMS Surface Roughness and Soil Textural 
Classification of Experimental Data Used 
to Develop Bare- Soil Algorithms 


Roughness Class 

Soil Textural Class 

RMS Height 
(cm) 

Number Of 
Data Points 

Smooth 

clay 

0.88 

13 


loam 

1,1 

16 


loam 

1.8 

17 


♦sandy loam 

1.32 

20 


♦silty '!;lay loam 

1.03 

20 


♦silty clay 

0.70 

20 

Medium Rough 

clay 

2.6 

13 


loam 

2.2 

17 


loam 

3.0 

19 

Rough 

clay 

4.3 

14 


loam 

4.1 

15 


* Data available at 10* and 20® In^^ldence angles only. 


89 





4.2.2 Vegeta tton-Covered Soli Algorithwis 


The form of the general algorithm for crops given by Equation 4.11 
assumes negligible two-way attentuatton of the soil backscatter component 
by the vegetation canopy. This assumption has been shown [14] to be 
valid for radar operating at frequencies less than 8 GHz and Incidence 
angles less than 20°. The breakdown of this assumption at higher angles 
of Incidence for Increasing canopy mass and transmission path length 
explains, In part, the decrease in least-square linear correlation co- 
efficient with Increasing Incidence angle as documented In Table 4.1. 

This decrease Is especially apparent for wheat and mllo which have pre- 
dominantly vertical canopy geometries. 

In lieu of extensive data sources for a° as a function of moisture 
for pasture and alfalfa at 4 to 5 GHz, the general algorithms given In 
Table 4.2 for these crops are derived from data for wheat. In addition, 
and for the same reason, the general algorithms for sand bars and mown 
pasture are assumed equivalent to that for smooth bare soil, 

The source data for the linear regression results given in Table 4.1 
are derived from measurements made over the full phenol ogic development 
of each crop and thus represent time-independent or mean canopy condi- 
tions. Also, the source data for each crop generally includes data from 
more than one experimental test field and more than one soil texture. 

All of the source data for the results shown in Table 4.1 were ob- 
tained from measurements made with no preferred radar look direction 
relative to crop row direction. As such, f(e) as given In Table 4.2 
Includes no effects due to periodic row geometry. 

From prior row-experiment measurements [43,44] for corn and soy- 
beans a simple look-direction modulation function is introduced to ac- 
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count for periodic soil structure underlying these crop canopies and 

mllo [13], It Is assumed t!v5it 5°/ from a look-direction parallel to crop 
row direction is equivalent to 5° from a vegetated surface without row 
patterns. The look- direction modulation function is defined by [13]: 

M(dB) • - aj/ (4.15) 

where 

» (j® In dB for radar look direction perpendicular to 
periodic ridge/furrow row structure, and 
■ of° in dB for radar look direction parallel to 
periodic ridge/furrow row structure. 

For corn and soybeans, whose mean values of M(dB) are 0.0 dB, 

0.5 dB, 1,0 dB, and 0,5 dB at incidence angles of 0°, 10°, 20°, and 
30° respectively, a third order polynomial fit gives M(dB) as a func- 
tion of Incidence angle. Thus, 

M{e) * 0.01670 + .0059^ - 0.167 X lO’^e^, dB (4.16) 

between 0° and 30° incidence angles. 

In Table 4.2 and for row crops with rows oriented perpendicular 
to radar look direction f(0)j^ “>s the sum of Equation 4.16 and f(9)// 
for the same crop aligned with row structure parallel to sensor look 
direction 

f(e)j_* f(e)// + M(0). dB (4.17) 

For simplicity, the above approach to row direction effects as- 
sumes: 

1) no dependence of M(e) on soil moisture. 
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2) < 3^11 Is equal to a° for a soil with no periodic row 
structure, and 

3) all three row crops have equivalent row structures 
and thus equivalent M(6). 

In addition, N(e) defined by Equation 4.16 applies only to dryland 
fanning characteristic of the simulation test site and characteristic 
of the source data. 

4.2.3 Target Classes with No Dependence on Soil Moisture 

Target classes from which radar backscatter at 4 to 5 GHz exhibit 
no dependence on soil moisture Include deciduous tree canopies, open 
water bodies, and cultural tar&4' , For the purposes of this simula- 
tion^ and in the absence of a comprehensive source data base, all cul- 
tural targets such as railroads, bridges, roads, and buildings are 
modeled as having a constant which will effectively saturate the 
simulated imagery. This approach Is consistent with Seasat Imagery 
for urban areas. 

Based upon data acquired by the University of Kansas [45], spring 
and sunmier deciduous tree canopies are considered to effectively atten- 
uate the backscatter component from the underlying soil at 4 to 5 GHz. 
Radar backscatter response to deciduous trees is therefore modeled 
as a function of cose for o° In natural units (m /m ) and 

a° trees » 10 log x cose), dB (4.18) 

Radar tackscatter from small water bodies is modeled after data 
acquired by the University of Kansas MAS 1-8 system In the spring of 1980. 
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These data were obtained at a frequency of 4.9 GHz for a pond 
roughly 1,5 acres in size with winds varying between 10 and 25 mph. 

The wind conditions are representative of conditions common to the 
simulated test area and produced waves of 1-to 2-cm in amplitude. 

A third-ordfir polynomial fit to data acquired at 0®, 10®, 20® and 
30® is assumed to be valid for wave conditions on the rivers and 
streams within the simulation. 

5° water * 22.82 - 5.126© + 0.237e^ - 3.973 x 1O’^0^, dB (4.19) 

4.2.4 Angular Dependence of 5° for Given Moisture Conditions 

For each target class, local incidence angle 9^, and soil moisture 
condition generated within the image simulations, a given scene element is 
characterized by a from Equation 4.11 and Table 4.2 as influenced by 
signal scintillation. Figures 4.2 and 4.3 illustrate the behavior of a° 
for select target classes as a function of e for MpQ equal to 10.0 and 
150.0 respectively. = 10 is the driest moisture condition permitted 

In the simulation and is roughly equivalent to the hygroscopic coefficient 
of a given soil, while Mp^ = 150 is the wettest moisture condition allowed 
and approximates a saturated soil condition. 

4.3 Simulation Data Base Construction 

The simulated data base corresponds to an 11 x 12 mile (17.7 km x 

19.3 km) rectangular area immediately to the east of Lawrence, Kansas. 

The Kansas River flows across the northern half of the simulated region. 

This region was mapped into a matrix of approximately 1,000 x 1,000 pixels 
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Mean Radar Backscattering Coefficient (dB) 



FIgurt 4 .2 Mean Baciacattering Coefficient of Several Target Classes os a 

Function of Angle of Incidence for Very Dry Soil Moisture Conditions, 

Mp^ ■■ 10.0% 
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Mean Radar Backscattering Coefficient (dB) 



* 150 « 0 % 
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where each matrix element represents a ground resolution of 18 x 18 
meters. A separate matrix was constructed for each of the following 
target parameters considered In the simulation: target category 
(Table 4.2), underlying soil texture, surface elevation, and Inci- 
dent rainfall from a simulated thunderstorm. 

4.3.1 Target Categories 

Target class data was extracted from color Infrared U-2 Imagery 
acquired in May, 1978. A positive color transparency of the 11 x 12 
mile area was rephotographed Into four quadrants as black and white 
negatives. Each quadrant could then be digitized Into a 512 x 512 
matrix on a CRT by use of a video camera. 

• Manual Interpretation of the U-2 imagery as augmented by low alti- 
tude USDA/SCS Imagery allowed assignment of a target class to each of 
the 512 X 512 pixel elements on each quadrant. A list of target cate- 
gories used to characterize land use within the data base is given in 
Table 4*2. Figure 4.4 shows the U-2 imagery used In interpretation 
and Figures 4.5 to 4.7 show the data base as classified for various 
target categories. 

Table 4.4 gives the percentage of the total 11 x 12 mile area as- 
signed to each target class in the data base. Pasture comprises the 
largest single crop class In the data base and is followed closely by 
corn, bare soil, deciduous trees, and soybeans. The percents of the 
total planted area for each crop class in the data base are compared 
in Table 4.5 to those reported for Douglas County, Kansas [46] as a 
whole in 1977 and 1978. 
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bridges, and buildings are shoMn 
in white. 



CRICIN.M 

BLACK AND WHITE PHOTOGKAPM 




■It 


TABLE 4.4 

Area Percent of Total Data Base 
Assigned to Each Target Class 


Target Class 

Percent of 
Total Area 

Roads 

3.76 

Railroads 

0.12 

River Bridges 

0.01 

City Structures 

0.85 

Rivers 

2.19) 

Lakes, ponds, 
impondments 

0.48) 

Smooth Bare Soil 

6.63\ 

Medium Rough Bare 
Soil 

14.30 

Rough Bare Soil 

2.75) 

Mown Pasture 

7.06 

Pasture 

15.93 

Alfalfa 

4.15 

Wheat 

6.65 

Sandbars 

0.35 

Deciduous Trees 

13.03 

Soybeans N/S Rows 
Soybeans E/W Rows 

5.46) 

5.62) 

Milo N/S Rows 

2.76) 

Milo E/W Rows 

2.27( 

Com N/S Rows 

8.32) 

Corn E/W Rows 

6,6l) 


N/S refers to crops planted with rows running north to south. 
E/W refers to crops planted with rows running east to west. 
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TABLE 4.5 

Comparison of Data Base Composition with 
Historical Records for Douglas County, Kansas 


Crop Class 

% of Total Area Planted 

Douglas County, 
Kansas [46] 

Simulation 

Data 

Base 

1977 

1978 

Com 

12.47 

16.05 

30.53 

Milo 

24.93 

16.32 

10.29 

Soybean 

12.47 

25.57 

22.66 

Alfalfa 

4.50 

7.07 

8.49 

Wheat 

23.55 

9.07 

13.60 

Hay (Mown Pasture) 

22.09 

25.93 

14.44 




For the two years shown, the crop>rn1x of the data base compares 
favorably to that reported for all of Douglas County with certain ex- 
ceptions. 309i of planted acreage In the data base Is corn, while only 
]ZX to 16% ii reported for the county In 1977 and 1978. Conversely, 
only 10% of the planted acreage In the data base Is mllo while 16%to 
25% is reported for the county as a whole. These discrepancies are 
expected since the Kansas River floodplain comprises a large portion 
of the data base. Much of the county's corn acreage Is concentrated 
within the floodplain, while most of the milo acreage Is In the upland 
areas typical of most of the county. 

4.3.2 Soil Textural Classification 

A digital matrix of soil textural class was constructed for all 
18 X 18 meter pixel elements In the data base. United States Depart- 
ment of Agriculture Soil Conservation Service county soil surveys of 
Douglas, Jefferson, Johnson, and Leavenworth Counties in Kansas were 
the data sources [47,48,49,50]. A photo-mosaic of SCS imagery inter- 
preted by SCS for mapping units of soil series and phase was re-inter- 
preted on an acetate overlay for A-horizon soil textural class. Each 
SCS soil phase was classified as having one of the 10 soil textural 
classes or complexes given in Table 4.6. The acetate overlay was then 
photo-reduced into quadrants and digitized on the CRT. Figure 4.8 
shows the digitized soil textural overlay for the simulation data base. 

Each soil textural class was assumed to have characteristic sand, 
silt, and clay components as given in Table 4.6. These textural com- 
ponents were then used to determine a characteristic l/3<*bar water 
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Son Texi 
Percent, an( 
Contents W1i 


Soil Textural Class 
Sand 

Loan\y sand 
Sandy loam 
Loam 

Silt loam 
Silty clay loam 
Silty clay 
Clay loam 

Complex (50^ loam and 
50% silty clay loam) 

Complex (50% silt loam 
and 50% silty clay 
loam) 


*1 /3-bar water content i 


TABLE 4.6 

j: 

i 

tural Classes, Their A'-ea • 

i Characteristic 1/3-^ ar Water I 

thin the Simulation Data Base I 


Textural Components 

1/3- Bar 
Water 
Content* 
FC 

Percent 
Area of 
Total 
Data Base 

% Sand 

% Silt 

% Clay 


5 

3 
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65 
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10 
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46 

.3375 

3.3 

33 
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33 
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25 
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26.5 

.2568 

0,7 
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61 

24 

.2723 
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Figure 4.8 Soil classes within the simulation test site. Pixel 
intensity is related to A-horizon soil texture from 
Table 4.6 where sand is darkest and complexes are 
brightest. The image «^epresents a sampling of every 
other column and row of the entire data base (25% sample). 
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content from Equation 4.12 for each soil textural class. The area 
percent of the total data base represented by each soil textural class 
Is also given In Table 4.6. 

4.3.3 Surface Elevation 

The distribution of surface elevation above mean sea level was 
encoded into the data base through procedures similar to those used 
for soil texture. In this case, the source data were United States 
Geological Survey 7.5 minute series topographic maps. An acetate 
overlay was made from a mosaic of the following quadrangles in Kansas: 
Lawrence East Quadrangle 
Baldwin City Quadrangle 
Eudora Quadrangle 
Edgerton Quadrangle 
DeSoto Quadrangle 
Gardener Quadrangle 

While these quadrangles have 10-foot contour intervals, it was generally 
necessary to construct the overlay with 20-foot contour intervals in the 
hill regions enveloping the river floodplain. In order to compute a 
realistic range to each piiiel, area of each pixel, and local angle of 
incidence, it was necessary to convolve the digital elevation data to 
produce local surface slope. The convolution procedure is treated in 
a later section. 

Figure 4,9 gives a shaded relief presentation of the final eleva- 
tion matrix after convolution. Intensity of each pixel element 
is related to elevation above mean sea leveV 
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Figure 4.9 Elevation above mean sea level. Every other column and 
row of the total data matrix are samples to produce 
this image. 


^•3*^ Data Base Registration 


The initial digital data base consisted of matrices for target 
class, soil texture, and surface elevation in four quadrants each. 

A set of 8 to 12 registration points had been preselected within each 
quadrant. An nth order polymomial fit to the registration points 
established a warping function for mapping all soil texture and sur- 
face elevation matrices into the coordinate system of the correspond- 
ing target class matrix. This process corrected for any lack of geo- 
metric equivalence in the raw data sources and any scaling differences 
arising from the digitization procedure. Upon completion of the warp- 
inon all four quadrants were then mapped into a common coordinate sys- 
tem, aligned, and overlap data were deleted from the composite data 
base. 

4.4 S imulation of Slope 

The elevation data base which was described previously consists of 
discrete heights in a matrix format. This matrix is operated upon by the 
simulation algorithms to produce slope information in the range and az- 
imuth dimensions, from which the local angle of incidence between the ra- 
dar antenna and the ground sample location is later computed. That is, 
for each and every entry in the ground-truth data base a facet orienta- 
tion is generated, with the terrain slopes between sample points modeled 
by linear interpolation of elevations in range and azimuth. 

If we define the angle of the slope in azimuth as p and the angle 
of the slope in the range directions as then the local angle of in- 
cidence 0 . is given as 
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(4.20) 


0 ■ cos“^ tan'j> sine + cose 
^ ‘ / tan^p + tan^\|) + 1 

assuming a side looking antenna. The angle of incidence e is defined as 
previously, in section 4.1.3, as the earth curvature corrected incidence 
angle. 

Taking several examples, we see that if and p are both small 
(level terrain cases), then e « e. If the ground is hilly, then fore- 
slopes (t|) positive) and backslopes (i|i negative) are present. For pos- 
itive and p small, the local angle of incidence can approach zero degrees 
and the for,‘.>lope backscatter can become quite large, For \J) negative 
and p small, the local angle of incidence can approach grazing (e. « 90®) 
and the backscattered signal has a small amplitude. The shadow case 
arises when a backslope has a value of |ij;l > e, and no signal is back- 
scattered from the shadowed terrain. 

In the situation of general terrain, the local slope treatment can 
be summarized by saying that a facet model slope in range and azimuth 
is calculated, based upon linear interpolation between elevation samples. 
Since power backscattered is computed on the basis of the resolution of 
the data base (as opposed to the sensor resolution, which may be greater 
than the 20-m data base resolution achieved herein) then the facets are 
smaller than the resolution cells of the sensor in general. 

A random model for tree heights was employed to achieve realistic 
conditions of natural growth. In the geometry calculations the random 
tree heights were incorporated such that the presence of taller trees 
could give rise to layover and shadow, which is commonly seen in very 
high resolution radar Imagery. That is, the random tree heights were 
calculated before processing of other geometric effects took place. 
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4.5 Simulation of Rainfall and Soil Moisture Conditions 


In orden to pMuce realistic simulations, four hypothetical soil 
moisture distributions were developed for the 0^5 cm soil layer. These 
distributions cover the full range of potential moisture conditions from 
saturation of the soil to a drought-like condition. With respect to 
time, the hypothetical moisture conditions cover a 35-day time span 
with simulated satellite overpasses on Days 4, 5, 15, and 35. 


4.5.1 Soil Moisture Conditions 

For each satellite overpass, the "actual" moisture of each 
18 X 18 meter pixel element In the data base Is established from a 
hypothetical rainfall and evaporative history. The objective Is to 
establish a set of moisture conditions for each radar simulation that 
Is a reasonable facsimile of common "real world" conditions^ In order 
to simplify the calculations, and because of a lack of adequate 
source data, several assumptions are made. First, the hydraulic con- 
ductivities of all soil textural classes In the data base are assumed 
to be equivalent. The bulk densities of all soils are assumed to be 

1.0 g/cm within the 0-5 cm soil layer. In addition, the evaporation 
rates from all soils and crop classes are assumed to be equal although 
still time dependent. In actuality of course, bulk density, evapor- 
ation rate, and hydraulic conductivity all vary as functions of many 
variables Including soil texture, crop class, and soil moisture. 

For each 18 x 18 meter pixel, soil moisture Me- is computed from 
the general expression: 


Hp- • Hp. + 100(M„ /FC„.,) 


(4.21) 
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where 


Mpc * 


M. 




FC 


voV 


% of 1/3 bar water content for the simulation 

the change In volumetric water content since the 
time of the (N-1) simulation, g/cn? 

estimated 1/3 bar water content, g/cm . 


Sincfi soil bulk density is assumed to be 1.0 g/cm^, then 


= 1.0 X FC 


( 4 . 22 ) 


for values of FC given in Table 4.6 for each soil texture. In addition, 

Mw Is dependent upon incident rainfall, surface and subsurface draln- 
''N 

age, and evaporation. 


'n ''(N-1) 


(4.23) 


where 


Mw = initial volumetric soil moisture 
''(N-1) 

3 

^rsin “ increase in My due to precipitation, g/cm 

Mevap " ^oss in My due to evaporative demand, g/cn? 

3 

Since soil bulk density is assumed equal to 1.0 g/cm , M^,.^ and M^„__ 

icn n cVap 

can be calculated from incident rainfall and evaporation rate by: 

Mrain “ (Precipitatic.i - Drainage) /Soil Depth (4.24) 


and 


Mgvap (Evaporation rate x elapsed time)/Soil Depth (4.25) 


where precipitation, drainage, and soil depth are in cm and evaporation 
rate is in cm per day. 
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4.5.2 Postulated Rainfall and Evaporative History 


During the first day of the 35-day time span, a steady and heavy 
rain Is assumed to have produced saturated so\H moisture conditions over 
the entire data base. Over the next three days, water In excess of 
field capacity (as estimated by the water retention of a given soil at 
1/3 bar) Is assumed to have drained from the upper 5 cm of soil through- 
out the data base. Thus a satellite overpass on Day 4 observes all 
soils with MpQ equal to 100^ of 1/3-bar water content, 

Less than a day after the first simulated radar overpass, a hypo- 
thetical convectional thunderstorm passes from west to east across the 
data base. The storm deposits a Gaussian rainfall distribution with 
maximum Incident rainfall of 2,5 cm along the center of the storm 
track and with a m1nimUi« approaching O.D cm along the northern and 
southern edges of the data base. Limited hydraulic conductivity of 
the soil causes all Incident rainfall in excess of 1.25 cm to drain 
laterally from the soil as surface runoff. Any incident rainfall less 
than 1,25 cm is assumed to percolate rapidly into the upper 5 cm of 
soil. Thus, immediately after passage of the thunderstorm, maximum 
Mrain 0.25 g/cm from Equation 4.24. A satellite overpass of the 
data base several hours after passage of the thunderstorm could be 
expected to observe a soil moisture distribution similar to the one 
defined above for Day 5 and used in the second radar simulation. 

During the following 10 days the upper 5 cm of soil dries due to 
evaporation at a rate of 0.1 cm of water per day. There is no additional 
rainfall. Thus Is a constant 0.2 g/cm over the data base from 

Equation 4,25. In addition, it is assumed that for any given pixel 
element Mep > 25.0% of 1/3-bar water content on Day 15. This assump- 
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tton reflects the capacity of most soils to replenish some portion of 
daytime evaporative loss during the night through capillary recharge 
and vapor flow. A lower limit of 2S% on surface moisture also functions 
to prevent unnaturally dry moisture conditions in the simulation which 
otherwise could approach or exceed "oven dryness" at 0^ Mp^ for sandy 
soils In the data base which have small volumetric water retentions at 
1/3 bar. Thus, a simulated radar overpass on Day 15 observes soil 
moisture conditions with a maximum of approximately 90X and a 
minimum defined at 25%, 

From the third simulation until the end of the 35-day period there 
is no additional precipitation, hence equals 0.0 g/cm . Since, 

at low surface-moisture conditions, evaporative rate is limited by the 
availability of near-surface water, the evaporation rate for this 20- 
day period is assumed to be reduced to 0.05 cm of water per day. As 
a result, in Equation 4.25 becomes 0.2 g/cm . In order to pre- 

vent MpQ of sandy soils from becoming zero, a lower limit of Mp^ - 10% 
is assumed to be valid on Day 35. Moisture conditions of 10% of 1/3- 
bar water retention approximates the hygroscopic coefficient of many 
soils and moistures less than this value are not readily attained under 
natural field conditions. Thus, a simulated radar overpass on Day 35 
observes moisture conditions which range between 32% and 10% of the 
1/3-bar water retention in the 0-5 cm layer. Such moisture levels 
are typical of those observed within the data base region during droughts 
in 1975 and 1976. 

The given moisture conditions for the four simulations as described 
above are tabulated in Table 4.7. Figures 4.10 to 4.13 show an image 
presentation of the moistures input into the simulations. Roads, build- 
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TABLE 4.7 

Values of Rainfall and Evaporation Used to Derive th 













# 
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igure 4.10 Input soil moisture for Moisture Figure 4.11 Input soil moisture for Moisture 

Condition 1, three days after a Condition 2, immediately after passage 

saturating rain. All soils are at ® thunderstorm. Surface soil moisture 

100% of field capacity. ranges fron field capacity to saturation. 
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black and whitc photograph 




Figure 4.12 Input soil moisture for Moisture Figure 4.13 Input soil moisture for Moisture 

Condition 3, ten days after the Condition 4, thirty days after the 

thunderstorm. thunderstorm, drought conditions. 



Ings, and water bodies are shown In black on these figures since soil 
moisture Is undefined for these categories and grey tone level Is pro- 
portional to MpQ. 

4.6 Generation of Radar Images - Summary of Computer 

Algorithms and Images 

The computer algorithms that produce the simulated radar Images, 
given the backscatter data, ground truth data base, and sensor parameters, 
are broken Into four major packages. They handle the following aspects 
of the process: (1) measuring the average value of the received power 

for each range gate, (2) Monte Carlo simulation of the fading, (3) aver- 
aging to trade spatial for radiometric resolution, and (4) scaling of the 
Irfiage data. We shall discuss these four steps briefly; the reader is 
further directed to ’[28] for greater detail. 

The simulation package mimics portions of the operation of a SAR for 
which digital processing is to be used. As such, a certain number of 
range gates are set up; in this case we wished to produce a ground range 
simulated image, so the gates are not of equal length across the image, 
(They monotonically increase in size with increasing incidence angle.) 
Next, the facet orientation calculations are performed for each point in 
the data base, and range to each facet is obtained. In order to deter- 
mine whether a location is in shadow, it is necessary to make neighbor- 
hood elevation comparisons, and this is done next. 

The "start” and "stop" ranges which mark the beginnings and endings 
of each range gate, having been determined previously, are available for 
approximating how much of each facet's return power will fall into a 
certain range bin. As the software exists at present, if the range to a 
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certain data base location equals the ending range of bin “K»'' then its 
power is equally divided between bins K and K f 1, and so forth for other 
proportions. The power referred to here is the average power determined 
by the radar equation [29] which accounts for antenna gains, diffrac- 
tion losses, scattering cross-section, and wavelength of illumination. 

Through the explicit computation of range bin stop and start ranges 
(or equivalently, times), the geometry/propagation factors of earth cur- 
vature, foreshortening, and layover aro automatically incorporated. An 
efficient binary search routine is actually used to determine when a 
facet's power will be received. 

There are certain measurable system parameters including transmitted 
power, antenna gain, wavelength, Doppler weightings, incidence angle, 
scale and bias factors, etc. in the SAR data processing that perhaps 
could be reJDOved in effect from the radar imagery that would allow one 
to use a soil moisture algorithm directly on the digital SAR image. Thus, 
the decision was made that the average received power would be calculated 
as 


c <,’(6^) A 


(4.26) 


where c incorporates all the factors excluded from the radar equation 
[29], and, furthermore, a reference range to the fourth power, to pre- 
vent P|^ from being very small. 

The above description has treated the problems of average power and 
geometry considerations, and concludes the discussion of the first of 
four packages of the simulation algorithms, as tailored to this particu- 
lar project. 


116 


The second of the four packages accounts for fading* which Is the 
major source of noise for small tine-bandwidth product sensors. The 
model derived in [33] is used to generate received power from average 
received power through Monte Carlo simulation as 

Pr * • y (4.27) 

where y "’s* y is described by a chi-square distribution with 

2N degrees of freedom. When a small number of independent samples is 

desired* (say, N < 5) Equation (4.27) provides a simple route for simulation 

» 

of the faded received power (the sum of the squares of the in-phase and 
quadrature voltage signals received). However, when large N images 
were needed, spatial domain smoothing to sacrifice resolution for looks 
was applied to a small N image. The spatial domain implementation of 
the continuous scanning mixed integration process [?9] was used for 
these types of images. For N = 12 images a Gaussian approximation to 
the multiplicative model of Equation (4.27) was used to produce values of 

K ' 

'’r * (4.28) 

where Z u(0,l), that is, Z is distributed as a Gaussian with mean zero 
and unity variance. 

The third of four packages of the simulation performs smoothing in 
either or both of the dimensions of the radar images. In several cases 
it was possible to use the same Image several times, at its finest pos- 
sible resolution and small N, and then at degraded resolutions and larger 
values of N. In all cases, the convolution weightings were obtained by 
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derivation of the filters whose structures are described and justified 
In [28]. j 

The fourth package of the simulation software simply took the values I 

of power produced by Packages 1» 2 and 3, and converted them to densities ! 

i 

through a logarithm operation, I.e., 0 » 10 log P^. This step was Incor- | 

porated simply to make the simulated Images compatible with the soil | 

moisture estimation algorithm. I 

A summary of the simulated Images Is given In Table 4.8. In all, i 

a set of 24 Images was produced for 3 resolutions x 2 angular ranges x 4 | 

y 

moisture conditions. Samples of these Images are given 1n Figures 4.14 to 
4.22 for selected conditions. Image graytone Gj Is scaled from 0 to 255 and j! 
Is proportional to received power In dB over a 50 dB dynamic range from 
-25 dB to -t25 dB. 

Gj • 5(10 log P^ 25) t 3 (4.29) 

Simulated Images for the 20-meter resolution case over the 7.5® to 9.3® 
angular range are shown In Figures 4.14 to 4.17 for moisture conditions 
of field capacity, saturation, drying, and drought, respectively. The drying 

i.; 

condition (Moisture Condition 3, 10 days after the hypothetical thunderstorm) I 

j* 

Is used as a case example of the effects of changing resolution and angular ( 
range In Figures 4.18 to 4.22. Because of matrix -size constraints on the I 

video display screen, all Images show a 25% sample of the total data base, j 

wherein every other column and row of the total digital matrix Is Imaged. ! 

As a result, certain linear features such as roads may appear to be discon- I 

h 

tinuous on the Images. | 

I 

5.0 EVALUATION OF RESlilTS | 

j 

The simulated satellite radar data was evaluated for its capacity to j 

correctly estimate the actual data-base moisture conditions given In | 

Section 4.5. Each of 24 simulations (2 angular ranges x 3 resolutions j 

X 4 moisture conditions) were interpreted by a generalized moisture 
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TABLE 4.8 

Suirmary of Simulated Radar Images 


N of 
looks 

Radar Pixel at 

Mid-Range Angle 

r. X r^, (m x m) 
a y 

r" 

Angular 
Range o 

12 

20 X 20 

7.5‘’-9.3^ 

12 

20 X 20 

ll.r-12.8° 

23 

93 X 100 

7.5‘’-9.3° 

27 

93 X 60 

ll.r=12.8‘’ 

2250 

1 ,000 X 1 ,000 

7.5®-9.3” 

2800 

900 X 600 

ll.r-12.8'^ 


All cases have simulated soil moisture conditions 
of; field capacity, saturation, drying, and drought. 
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Figure 4.14 Radar simulation with 20-meter Figure 4.15 Radar simulation with 20-meter resolution and 

resolution and 7.5“ angle of 7.5“ angle of incidence for Hoisture Condition 

incidence for Moisture Condition 1, imnediately after a thunderstorm, 

all soils ai, field capacity. 
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Figure 4.18 Radar simulation with 100-meter Figure 4.19 Radar simulation with 1-kilometer 

resolution and 7.5° angle of incidence resolution and 7.5° angle of incidence 

for Moisture Condition 3, 10 days for Moisture Condition 3, 10 days after 

ai ter the thunderstorm. the thunderstorm. 
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igure 4.20 Radar simulation with 20-meter ^solution and 4 2i p^Jdar simulation with 100-meter resolution and 

11.1 angle of incidence for Moisture Condition 3, 11.1® angle 0 ^ incidence for Moisture Condition 

10 days after the thunderstorm. 10 days after the thunderstorm. 
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angle of incidence for Moisture 
3, 10 days after the thunderstorm. 


algorithm A(e) to estimate for each resolution cell. Estimated 
moisture MpQ was then correlated with actual moisture to establish 
estimate accuracy for a given simulation, This was accomplished using 
a digital processing scenario for the simulated radar data, although 
Image processing of the quantized radar data should produce similar 
results. 

Evaluation of sensor performance in terms of moisture estimate 
accuracy proceeded at several levels dependent upon the assumed avail- 
ability of ancillary data, Ancillary data includes: 

1. known location of water bodies, buildings, roads, and bridges; 

2. a pHori information about whether an agricultural scene is 
bare soil or covered with a crop canopy. 

Knowledge of position and extent of water bodies, buildings and roads 
could be the product of a manually or machine-interpreted image from 
one of the simulated radars, some other image source, or digital map 
data base. Agricultural crop-cover information may be obtained through 
the use of local crop calendars or may be obtained from crop discrimination 
from Landsat and/or a satellite radar operating at 14 GHz with 50° 
incidence angle. 

For the simple case where there is no ancillary data concerning 
the imaged area, the simulated backscatter data is interpreted using 
a blind classifier to estimate moisture A(e). In this case, all simu- 
lated image elements are interpreted by the generalized algorithm 
as containing valid surface soil moisture information. 

In the second case, where the positions of water bodies, buildings, 
and roads are assumed to be known, the pixel elements representing 
these l^atures are filtered from the simulated backscatter data and therefore 
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not assigned a soil moisture estimate Mp^. In this approach, the 
evaluation of estimate-accuracy is less biased by the misinterpretation 
of cultural features and water bodies as agricultural scenes, since 
the generalized moisture algorithm pertains only to agricultural scenes 
and soil moisture is undefined for water bodies and hard targets. 

The third case, in its simplest form, permits the division of agri- 
cultural-scene data into that representing bare soil and that represent- 
ing any crop canopy. Of course, if complete crop discrimination were 
assumed as auxiliary data then the return from each identifiable crop 
cover category could be treated independently. The simulation evalua- 
tion procedure treated herein covers only the first two cases and the 
simple form of the third where two generalized algorithms are used to 
interpret the backscatter data, one algorithm Q(^) for generalized 
bare soil conditions and one algorithm C(e) for generalized crop canopy 
conditions. 

5.1 Generalized Moisture Interpretation Algorithms 

The generalized moisture algorithm A( 6 ) was developed for an amal- 
gamation of all agricultural-scene data from sources as listed in 
Table 4.1. The available data were at frequencies between 4.25 GHz 
and 4.9 GHz and incidence angles between 0° and 30®. In addition, 
these same data sources were used to develop general moisture interpre- 
tation algorithms for all bare soil classes 6 ( 0 ) as considered separ- 
ately from all crop canopy classes 0 ( 0 ). 

The amalgamated data at each angle (0®, 10®, 20°, and 30*) was 
submitted to least-squares linear regression to establish the linear 


126 


dependence of cr® on The regression results are shown in 
Table 5.1 for all classes canbined, bare soil classes only, and 
vegetation canopy classes only. 

Scattergrams of a° as a function of «FC for all agricultrual 
scene data combined are shown in Figures 5.1 and 5.2 for 10" and 20® 
incidence angles respectively. Linear correlation coefficient 9 
is observed to decrease from 0.88 at 10® to 0.65 at 20® and Is felt 
to be caused by the increasing impact of canopy attenuation and 
canopy geometry at angles of 20® and greater. 

Assuming normal distributions of o° and in Figure 5.1, the 
90^ confidence Interval about is computed to vary from +/- 0«20 dB 
at MpQ = 75 to a maximum of +/- 0.53 dB at Mp^ = 0 or 150. This means 
that 90% of the time Mp^ should be within +/- 1*5% of the measured Mp^ 
when Mpj, = 75% and ftp^. should be within +/- 3.9% of the measured Mp^. 
when MpQ equals 150% or 0%. In a similar fashion, calculation of the 
90% confidence interval about the linear regression at 20® 

(Figure 5.2) yields a .9 probability that measured Mp^ is within ±3.2% of 
MpQ for Mp^ = 75% and within +/- 8.7% of MpQ for the moisture extremes 
where Mp^ = 0% or 150%. 

Scattergrams of all bare soil cases are shown in Figures 5.3 and 
5. A for 10® and 20® respectively. These plots include 181 data points 
from 11 test fields where roughness ranges from RMS heights of 0.9 cm 
to 4.3 cm and soil texture ranges from sandy loam to clay (Table 4.3). 
Table 5.1 shows that for bare soil cases considered alone, linear cor- 
relation coefficient maximizes at 10® and decreases markedly at nadir 


TABLE 5.1 


Results of Least-Square Linear Regressions 
Used to Establish the General Interpretation 
Algorithms for: (1) All Classes, (2) Bare Soil Classes, 
and (3) Crop Canopy Classes 


Target Classes 

Sample 

Size 

Inci- 

dence 

Angle, 

0 

Regression 

Coefficients 

Linear 

Corre- 

lation 

Coeffi- 

cient 

Inter- 

cept, 

dB 

Slope 

dB/1,0^ 

"fc 

All Classes Combined 

268 

0 

- 9.67 

0.162 

.651 


324 

10 

-14.34 

0.133 

.883 


324 

20 

-14.80 

0.078 

.653 


268 

30 

-16.01 

0.068 

.505 

All Bare Soil Classes 

124 

0 

-10.92 

0.170 

.486 


181 

10 

-15.96 

0.148 

.849 


181 

20 

-17.86 

0.114 

.700 


124 

30 

-20.61 

0.121 

.588 

All Vegetation Canopy 

144 

0 

- 9.38 

0.165 

.855 

Classes 

143 

10 

-13.84 

0.133 

.921 


143 


-13.03 

0.048 

.590 


144 

30 

- 

-14.33 


.525 


128 




ORIGINAL PAGE IS 
OF POOR QUALITY 



^ 00 

CSJ •— I 


CM 


nO 


(8P)oO 


NO 

I 


Csl 

I— I 

I 


00 

rH 

I 


I 


129 


Percent or 1/5 Bar Water Content in the 0-5 cm Layer 

Figure 5.1 Backscatter coefficient as a function of Mpc in the 0-5 cm soil layer 
for 324 observations of bare and vegetated soils at 4.5 GHz, HH polari- 
zation, and 10'" incidence angle. 



078Mcr-14.80 



Percent of 1/3 Bar Water Content in the 0-5 cm Layer 

Backscattering coefficient as a function of Mpc in the 0-5 cm soil 
layer for 324 observations of bare and vegetated soils at 4.o faHz, 
HH polarization and 20® angle of incidence. 
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Percent of 1/3 Bar Water Content in the 0-5 cm Layer 

Figure 5.3 Backscattering coefficient as a function of Mpc in the 0-5 cm soil 
layer for 181 observations of bare soil at 4.5 6Hz, HH polarization 
and 10° angle of incidence. 
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Figure 5.4 Backscattering coefficient as a function of Mpc in the 0-5 cm soil 
layer for 181 observations of bare soil at4.5 GHz, HH polarization 
and 20® angle of incidence. 



and as incidence angle approaches 30®. This behavior is caused by the 
effects of surface roughness near nadir and at high incidence angles. 
Assuning normal distributions, the 90?l{ confidence interval for the 
bare soil data plotted in Figures 5.3 and 5.4 gives: 


"fc 

90 % Confidence Interval dB 

10® 

20® 

25 

+/- .74 

+/- .94 

75 

+/- .31 

+/- .39 

100 

+/- .33 

.42 

150 

+/- .78 

+/-1.0 


Scatterplots of 143 data points from fields of corn, soybean, milo, 
and wheat are shown in Figures 5.5 and 5.6 for 10® and 20®, respectively. 

For each crop these data cover the full range of crop growth stages and 
also represent different soil textures and roughness conditions (Table 4.3). 
The linear correlation coefficient decreases from 0.92 at 10° to 0.59 
at 20® because of increasing interclass variance at higher angles due to 
surface roughness and canopy attenuation differences and also because of 
higher intraclass variance within the cases of wheat and milo. Ninety percent 
confidence intervals for crop canopy data alone in Figures 5.5 and 5.6 
are similar to those given above as computed for bare soil at 10° and 20°. 

For each of the three sets of agricultural scene classes (all 
classes, bare soil only, crop canopy only), the linear regression coefficients 








Vegetation Data 
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Figure 5.6 Backscattering coefficient as a function of MpQ in the 0-5 cm layer 
for 143 observations of various crops at 4.5 GHz, HH polarization, 
and 20° angle of incidence. 



given In Table 5.1 were fitted with a third-order polynomial 
function of incidence angle. Thus, intercept and slope for each gen- 
eral scene class are expressed in terms of f(e) and g(e) respectively 
as given in Table 5.2 and shown in Figures 5.7 and 5.8. These general 
equations of the form: 


a° * f(e) + g( 0 )M| 


FC 


(5.1) 


are then inverted to estimate as a function of simulated o°. 

Mpc * (a° - f(0))/g(9) (5.2) 

Thus, for the general case where all agricultural scene data is con- 
sidered simultaneously, 

A(e) * (a° + 9.67 + 0.840 - 4.59 x lO'^O^ + 8.27 X 10"%^^/ 
(0.161 + 9.38 X 10“% - 4.97 X 10“%^ + 1.21 x 

10“%% (5.3) 


Where 

A 

A(o) * estimated moisture MpQ as a function of incidence 
angle 0 relative to the satellite radar. 

A(e) is used to estimate Mp^ from simulated satellite radar data in the 
case where no ancillary data is available and also in the case where the 
locations of roads, buildings, and water bodies are assumed to be known. 

n 

In the case where ancillary data allows difffi(^ntial treatment of 
agricultural scene data into that from croj^ canopies and that from bare 
soil surfaces, moisture is esttmated from B(e) and C(e), respectively. 
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TABLE 5.2 

Third^Order Polynomial Fits of Linear Regression Coefficients 
As a Function of Incidence Angle Between 0® and 30® for 
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.Vegetation 

^ 

All Classes 


10 15 
Angle of Incidence (Degrees) 


Intercept of general i 2 cd algorithms as a function of angle of 
Incidence between 0® and 30®. f(e) Is shown for (a) all agri- 
cultural scene data combined, (b) bare soil data only, and (c) 
crop canopy data only. 




Sensitivity d6/1.0%M 
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Figure 5.8 Slope of generalized algorithms as a function of angle of 

incidence between 0° and 30°. g(e) is shown for (a) all agri- 
cultural scene data combined, (b) bare soil data only, and 
(c) crop canopy data only. 


B(0) « (o® 10.92 + 8.37 x 1O“'*0 - 4.06 x 1O"^0^ 

+ 7.84 X 1O‘^0^)/(O.17 + 6.02 x 10’% 

- 3.76 X 10‘V + 1.00 X 10“V) (5.4) 

C(0) * (a° 9.38 + 9.57 X 1O’^0 - 6.34 x lO’V 

+ 1.23 X 1O*^0^)/(O.165 + 4.0 X 1O’^0 

- 9.47 X lO'S^ + 2.27 X lO’V) (5.5) 


where 

B(0) ■ estimated moisture MpQ In the 0-5 cm layer 
of bare soil 

C(e) * estimated moisture Mp^ in the 0-5 cm layer 
beneath crop canopies. 

No attempt is made to account for the effects of crop type or row 
direction effects in C(e). 


5.2 Interpretation of Simulated Radar Images 

The general interpretation procedure is outlined in Figure 5.9. 
For a given radar simulation (angle swath, resolution, and moisture 
condition) the range Information associated with each o° value is 
used to calculate an effective incidence angle 0 relative to the radar. 
This computation assumes spherical earth geometry, constant orbital 
altitude relative to mean sea level, and a constant mean elevation of 

tv 

the data base. Estimated soil moisture in the 0-5 cm layer MpQ was 
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INTERPRETATION PROCEDURE 



Figure 5.9 General interpretation procedure for evaluating soil 
moisture estimate accuracy. 

141 



calculated from the general interpretation algorithm (Equation 5.3) 
given a° and 0 for each simulated radar image pixel element. Estimated 
moisture images were produced for each simulation by converting Mp^ 
to gray tone level by: 

S * 3/2 Mpj. + 16 (5.6) 

where 

Gpi " gray tone level between 0 and 255. 

Equation 5.6 also describes the scaling used to produce the "actual" 
moisture images in Section 4.5. 

Examples of estimated moisture images at the 7.5° to 9.3° angular 
range and for each resolution and moisture condition are shown in Fig- 
ures"5.10 to 5.15. All images have equivalent hue, intensity, and 
saturation scaling. The estimated soil moisture ”fc is shown in 
Figures 5.10 to 5.15 as derived from Equation 5.3 assumes no natural 
limits on estimated moisture, whereas real-world conditions limit 
"actual" soil moisture to 056 - Mp^ - 150^^. In addition. Equation 
5.3 represents a blind moisture classifier since it requires no 
a knowledge of local slope, land-use category, crop canopy, 

crop row orientation, or surface roughness. 

5.3 Evaluation of Moisture Estimate Accuracy Using Blind Classifier 

Precise evaluation of soil moisture estimate accuracy was complicated 
by the geometric relief displacements inherent in the radar imaging pro- 
cess. While the relatively small vertical relief and slope of the data 
base produced insignificant layover and shadowing in the resultant images. 
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Figure 5.14 Predicted moisture map estimated by Figure 5.15 Predicted moisture map estimated by 

general algorithm from radar general algorithm from radar 

simulation with 100-meter resolution simulation with 1-kilometer resolution 

and 7.5® angle of incidence for and 7.5® angle of incidence for 

Moisture Condition 3. Moisture Condition 3. 


foreshortening produced a geometric displacement of upland surfaces by as 
much as 100 meters relative to the surface defined by the river flood- 
plain. 

Quantitative evaluation of mol sure estimate accuracy by machine 
assumes that the radar Image can be accurately mapped back Into the 
coordinate system defined by the original data base. Because of the 
complexity of Image rectification to account for range creep caused 
by foreshortening In the radar images a dual appraoch was followed to 
establish estimate accuracy. In the first and 'Vorst-case" approach, 
moisture estimate error is examined for the entire data base ar.d Is subject 
to errors associated with Imaging geometry problems and the effects of 
unknown local slope. In the second approach, moisture estimate accuracy 
will be considered only for the relatively flat region defined by the 
river floodplain where elevation of the data base is less than or equal to 
820 feet. Both analytic approaches utilized a common set of moisture 
estimate error maps. 

Estimate error maps were produced by the following procedure. A 
set of five control points were identified on the simulated radar imagery 
for recognizable features on the floodplain. These same features were 
located on the data base. The mean distance between the control points 
as identified on the radar Imagery and those on the data base was used 
to translate the coordinate system of the estimated moisture image 
without rotation or warping. The resultant estimated moisture matrix was 
then compared to the “actual" moisture matrix from the data base with a 
mean registration error of +/- 20 meters on the floodplain and +/- 100 meters 
on the upland surfaces due to uncorrected range creep. Thus for a given 
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resolution and moisture condition, the estimated and ''actual” Mp^ 
could be compared on a pixel by pixel basis relative to the 18 x 18 meter 
resolution cells of the data base. 

Figures 5.16 to 5.19 show the difference between estimated and 
actual moisture for each resolution and moisture condition. 

Op * Mpp - Mpp + 128 (5.7) 

^1j '^^Ij 

where 

Gp K Image graytone level between 0 and 255 

Mpp * estimated percent of field capacity In the 0-6 cm layer 

Mpp •* "actual” percent of field capacity In the 0-5 cm layer 

1j “ Image position of a given pixel elerpcnt 

Thus, for a given image « 128 (medium gray) represents zero differ- 
ence between estimated and actual moisture, bright pixel elements with 
Gp » 128 represent a large overestimate of moisture, and dark pixel 
elements with Gp « 128 correspond to large underestimates of "actual” 
moisture. 

For any of the three simulated resolutions, some general observa- 
tions can be made regarding estimate error as a function of certain data- 
base characteristics. 

1. Pixels with Zero Moisture in the Data Base . By definition, 
all target classes in the data base for which "actual" moisture jlis un- 

i - 

defined and arbitrarily set to zero will result in large moisture 
estimate errors. Bridges, buildings, and railroads with a constant 
inean a° of 10.0 will always produce a moisture estimate » zero. For 
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Figure 5.18 So’l moisture estimate error, predicted Figure 5.19 Soil moisture estimate error, predicted 

minus input radar minus input radar 

simulation with 1-kilometer resolution simulation with 100-meter resolution 

and 7.5° angle of incidence for and 11.1° angle of incidence for 

Moisture Condition 3. Moisture Condition 3. 


water bodies and roads, "actual” moisture is undefined and set equal 
to zero, while estimated moisture for these target classes is always 
greater than zero for e near nadir. Thus, such targets are generally 
bright on the error maps (Figures 5.16 to 5.19)- These targets com- 
prise 7.45^ of the data base. 

2' Tree Canopies . Tree canopies are assumed to completely atten- 
uate the backscatter contribution from the underlying soil at the 
simulated frequency. Thus, estimated moisture is generally far less 
than "actual" moisture except for the very dry conditions simulated 
for Moisture Condition 4 (30 days after the thunderstorm). Deciduous 
trees comprise 13.0^ of the data base. 

3. Local Slope Effects . Since Equation 5.3 is a blind classifier 
of the radar image, a predictable error componeni is irit'^oduced into the 
error maps because of the high sensitivity of,o° to 0^ near nadir. 

In general, moisture is overestimated for west-facing local slopes 
(those toward the satellite) while moisture is underestimated for east- 
facing local slopes (those away from the satellite). As would be ex- 
pected, these effects are most noticeable for the dissected upland 
areas to the north and south of the river floodplain and for 20-meter 
resolution. Coarser resolution tends to average many of these errors 
associated with the effects of local slope. 

Range Creep . The translation of features on the radar imagery 
relative to the data base as a function of elevation by foreshortening 
produces a double error component on the error maps. Since radar 
image to data base rectification, was performed to minimize position 
errors on the floodplain, this error component is most noticeable in 
the dissected upland region in the lower right side of the error maps. 


150 


Moisture estimate errors in this region resemble "ghost images," 
especially for linear features such as roads. The "ghost" errors are 
approximately equal in magnitude but opposite in sign and are separated 
by one to five pixel elements. These double estimate errors are arti- 
facts of the comparison methodology and relate to the position of a 
feature such as a road on the data base and its offset location on 
the interpreted radar imagery. 

Figures 5.20 to 5.22 are the result of color slicing the graytone 
scale on the error maps shown in Figures 5.16 to 5.18 at the leveh 
given in Table 5.3. The error map pixels encoded in yellow corres- 
pond to an estimate error of +/- 20% of the "actual" input percent 
of l/3'-bar water content. This level of error in percent of field 
capacity represents an equivalent uncertainty in gravimetric moisture 
of +/- 1,8%, +/- 4,2%, and V- 5.1% for loamy sand, loam, and clay 
loam respectively from Table 5.3 

The moisture estimate accuracy as presented in Figures 5.20 to 5.23 
can also be related to initial estimates of accurai.y level required 
by agronomists and hydrologists. Tables 5,4 and 5.5 present the 
accuracy requirements identified for the soil moisture information 
user community [51]. The approximate relationship between accuracy 
levels in Tables 5.4 and 5.5 and estimate error in percent of field 
capacity is given in Table 5.6. 

5.4 Analysis of Moisture Estimate Accuracy Over the Total Data Base 

The relationships between absolute estimate error and cumulative 
percent of the 800,000 pixels in the total base are presented in 
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Figure 5.20 Soil moisture estimate error Figure 5.21 Soil moisturo estimate error 

level-sliced at accuracy levels level-sliced at accuracy 

3 and 4 from radar simulation with 3 and 4 from radar sinKilatlon with 

20-meter resolution and 7.5“ angle 100-meter resolution aid 7.5 angle 

of incidence for Moisture Condition 3. incidence for Mc’stjre Condition 
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Figure 5.22 Soil moisture estimate error 

level-sliced at accuracy levels 
3 and 4 from radar simulation with 
1-kilometer resolution and 7.5“ angle 
of incidence for Moisture Condition 3 
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TABLE 5.3 

Interpretation of Soil Moisture Estimate Error Maps 
E = Mpc - «FC 

where E = error in % of 1/3- bar water content 

r\ 

Mpp = % of 1/3-bar water content as estimated from radar image 
using Equation 5.3 

rlpQ = "actual" % of 1/3-bar water content from Equation 4.21 
A. Error Map Color Codes 


Color 

E, percent 

E max 

E min 

Blue 

— 00 

- 30 

Green 

- 30 

- 20 

Yellow 

- 20 

+ 20 

Orange 

+ 20 

+ 30 

Red 

+ 30 

+ 00 


B. Relationship of E to Uncertainty in Gravimetric Moisture Content 


|E| 

Gravimetric Moisture 

loamy sand 

loam 

clay loam 

30% 

+/- 2.69% 

+/- 6.33% 

+/- 7.63% 

20% 

+/- 1.80% 

+/- 4.22% 

+/- 5.09% 

10% 

+/- 0.90% 

+/- 2.11% 

+/- 2.54% 
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Soil Moisture Information and Data Requirements 
at Different Crop Production Stages (from [51]) 
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General accuracy of High, Medium, or Low 
Gradation between Accuracy Levels 1 and 5 
+/- 4% accuracy by value measurement 









TABLE 5.5 

Soil Moisture Information and Data Requirements in Hydrology (from [51]) 
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General accuracy of High, Medium, or Low 
Gradation between Levels 1 and 5 
±2% accuracy by volume raeasuranent 



TABLE 5.6 

Approximate Relationship Between Uncertainty in 
Percent of Field Capacity and User Require- 
ments of Moisture Estimate Accuracy 


User- Defined 
Accuracy 
Level [51] 

Uncertainty in 
% of Field 
Capacity 

1 

+/- 50% 

2 

+ 

1 

o 

3 

+/- 30% 

4 

+/- 20% 

5 

o 

1 

+ 



Figures 5.23 to 5.26 for Moisture Conditions 1 to 4 respectively for the 
7.5® to 9.3® angular swath. Similar plots for the 11.1° to 12.8® angular 
swath show only slightly lower estimate accuracy for each simulated mois- 
ture condition and only Moisture Condition 3 for drying conditions Is 
shown In Figure 5.27. These results represent a "worst case" evaluation 
of moisture estimation accuracy since no adjustment is made to account 
for errors caused by geometric registration problems, the effects of 
local slope, or the Inclusion of data base categories where moisture Is 
undefined. 

For Moisture Condition 1 where the soil Is at 1003^ of field capacity, 

A 

a +/- 20 % estimate accuracy Is achieved (estimated moisture M^q ranges 
between 80^ and 120%) over 55.6%, 65.0%, and 86.6% of the data base for 
resolutions of 20 m, 100 m, and 1 km respectively. This trend toward 
Increasing estimate accuracy at coarser resolutions also Is apparent In 
Table 5.7 and 5.8 for simulation Moisture Conditions 2 and 3 (Immediately 
after the thunderstorm and 10 days later). This result Is apparently 
due to several considerations: 

1. From Equations 4.21 to 4.25, "actual" moisture does not vary 
dramatically between adjacent 18 x 18 meter data base pixels, except 
at the boundaries of soil types, 

2. the coarser resolutions average the local effects of slope, 
canopy cover type, row direction and surface roughness, and 

3. the effects of water bodies, cultural targets, and forested 
areas are averaged over much larger areas. 

For Moisture Condition 4, drought conditions, the lower accuracy of 
the 1-km resolution compared to the finer resolutions is caused by the 
effects of cultural targets. Categories such as buildings, bridges. 





Figure 5.23. Soil Moisture Estimate Accuracy for the Entire Simulation 
Data Base with Moisture Condition All Soil at Field 
Capacity^ and at Angle of Incidence Between 7,5 and 9.3 . 
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100 meters 


Moisture Condition :02 
Angie of incidence 7.5® to 9.3® 
800,000 Pixels m Data Base 


Absolute Difference Between Actual and Predicted Mp< 

Figure 5.24, Soil Moisture Estimate Accuracy for the Entire Simulation 
Data Base with Moisture Condition 2 , Saturated Soil, 
and at Angles of Incidence Between 7.5® and 9.3®, 
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%of Database with Less Error 





Figure 5,25. Soil Moisture Estimote Accuracy for the Entire Simulation 
Doto Base With Moisture Condition 3, 10 Days After ^ 
Thunderstorm/ and at Angies of Incidence Between 7,5 

and 9,3®. 


161 



% of Database with Less Error 


ORIGINAL PAGE IS 
PO 


100 meters 


!r 


1 km 


20 meters 


Moisture Condition*. 04 
Angle of I ncidence 7. 5® to 9. 3® 
800, 000 Pixels in Data Base 


0 20 40 60 80 100 

Absolute Difference Between Actual and Predicated M 

Figure 5.26. Soil Moisture Estimate Accuracy for the Entire Simulation 
Data Bose with Moisture Condition 4, Drought, and at 
Angles of Incidence Between 7,5® and 9.3°. 


162 





ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 5.27 Soil moisture estimate accuracy for the entire simulation 
data base with Moisture Condition 3, 10 days after thunder 
storm, and at angles of incidence between 11.1® and 12.8®. 
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TABLE 5.7 

Cumulative Percent of Total Data Base Within Certain 
Absolute Estimate Error Limits for C-Band 
Simulation Results at e « 7.5® - 9.3® 


Approx. User Accuracy Level 
[51]-. 

1 

2 

3 

B 

5 

+/- % of field capacity 

50 

40 

30 

20 

10 

Moisture 

Condition 

+ 

Resolution 






1 

20 X 20 m 

79.2 

j 

74.5 

67.2 

55.6 

34.8 


93 X 100 m 

86.4 

82.7 

76.9 

65.0 

42.7 


1 X 1 km 

92.5 

92.1 

90.7 

86.6 

70.1 

2 

20 X 20 m 

76.3 

70.3 

61.3 

47.6 

27.6 


93 X 100 m 

84.2 

79.8 

71.7 

58.3 

34.8 


1 X 1 km 

92.4 

92.1 

90.8 

85.9 

61.3 

3 

20 X 20 m 

87.5 

81.3 

f 

72.5 

58.9 

35.9 


93 X 100 m 

92.0 

i 87.7 

81.0 

68.6 

44.2 


1 X 1 km 

92.9 

91.0 

i 

87.4 

80.4 

58.6 

4 

20 X 20 m 

93.8 

89.6 

80.9 

62.3 

1 34.0 


93 X 100 m 

93.9 

91.5 

86.5 

72.5 

, 42.1 


1 X 1 km 

87.3 

83.5 

79.3 

72.6 

1 50.3 
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TABLE 5,8 

Cumulative Percent of Total Data Base 
Within certain Absolute Estimate 
Error Limits for C-Band Simulation 
Results at e ■ 11.1®-! 2. 8® 


1 

1 






Approx. User Accuracy Level -► 

1 

2 

3 

4 

5 

+/- % of field capacity ^ 

50 

40 

30 

20 

10 

Moisture 
Condition + 

Resolution ' 1 ' 






1 

20 X 20 m 

79.0 

74.2 

67.2 

54.9 

33.0 


93 X 60 m 

83.7 

79.3 

72.7 

60.8 

37.3 


900 X 600 m 

91.0 

89,8 

87.1 

80.6 

62.8 

2 

20 X 20 m 

74.8 

70.1 

62.2 

49.3 

29.2 


93 X 60 m 

80.6 

76.1 

68.1 

54.7 

33.8 


900 X 600 m 

91.8 

91.0 

88.6 

81.8 

55.6 

3 

20 X 20 m 

89.0 

83.1 

73.6 

57.8 

33.1 


93 X 60 m 

92.0 

87.6 

79.3 

64.6 

38.9 


900 X 600 m 

91.1 

88.7 

84.5 

75.7 

53.8 

4 

20 X 20 m 

92.4 

86.8 

76.0 

57.8 

32.7 


93 X 60 m 

93.6 

90.2 

82.8 

67.3 

37.6 


900 X 600 m 

83.4 

81.2 

79.0 

74.2 

49.5 















and railroads have a constant of 10 dB and for very dry moisture con- 
ditions the backscatter from agricultural terrain adjacent to such tar- 
gets is typically 20 to 25 dB less. Thus during extremely dry conditions, 
the large spatial averaging inherent in the 1 km resolution will cause 
a significant overestimation of the moisture present in agricultural scenes 
adjacent to cultural targets characterized by a very large a°. As a con- 
sequence, it is expected that the spatial density of such targets within 
a given area will effectively determine the upper limit of desireable 
resolution for accurately sensing very dry soil moisture conditions. Thus, 
for agronomic regions similar to the data base, resolutions on the order of 
- 1 km would be adequate for dry conditions, while for areas such as the 
high plains (where the density of cultural targets is lower than that 
simulated) resolution t 1 km might prove adequate. On the other hand, 
accurate sensing of very dry moisture conditions in regions with dense 
distributions of hard cultural targets such as the northeastern United 
States and northern Europe would require a resolution less than 1 km. 

The estimate accuracy of the 20 m resolution is seriously affected by 
local slope and crop canopy cover. The effects of local slope are most pro- 
nounced for high moistures (Moisture Conditions 1 and 2) since the sen- 
sitivity of a° to local incidence angle is greatest for high soil mois- 
ture. This is reflected in the generally lower estimate accuracy of the 
20 meter resolution for Moisture Conditions 1 and 2 than for Moisture 
Conditions 3 and 4 where the angular dependence of is less pronounced. 

5.5 Analysis of Moisture Estimate Accuracy Within the Floodplain 

In order to minimize the analytic bias introduced by uncorrected 
foreshortening of the predicted moisture maps, moisture estimate accuracy 
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Is presented in Tables 5.9 and 5.10 for each resolution and moisture 
condition for only those pixels on the relatively flat river floodplain. 
Approximately 183,000 data-base pixels con^rise this region or about 
23 % of the data base. Geometric registration of the predicted moisture 
maps to the coordinate grid of the data base was optimized for this 
region. The floodplain was arbitrarily defined by all pixels in the 
data base with elevation less than or equal to 820 feet. In addition, 
the values given in Tables 5.9 and 5.10 exclude pixels classified in the 
data base as water bodies, buildings, railroads, bridges, roads, trees, 
and sandbars, thereby minimizing that component of net estimate error 
related to pixels with undefined "actual" soil moisture. 

While the values in Tables 5.9 and 5.10 are significantly higher 
than those shov»!n in Tables 5.7 and 5.8 for the entire data base, they 
reflect the same general trends with respect to moisture condition and 
resolution. The 1-km resolution results in the highest estimate accuracy 
at all accuracy levels except for the very dry moisture condition pre- 
sent 30 days after the thunderstorm. Also, estimate accuracy is gen- 
erally lower at all resolutions for either the extremely wet or extremely 
dry moisture conditions (Moisture Conditions 2 and 4 respectively) than 
for the more typical intermediate moisture conditions. This is expected 
since the effects of surface roughness and crop canopy cover cause a 
divergence of 5° between the target classes at the moisture extremes. 

. The cumulative percent of the 183,000 pixels comprising the agri- 
cultural floodplain is plotted versus maximum absolute estimate error 
in Figures 5.23 to 5.31 for Moisture Conditions 1 to 4, respectively, 
for the 7.5® to 9.3® angular swath and in Figure 5.32 for Moisture 
Condition 3 and the 11.1® to 12.8® angular swath. At Accuracy Level 4, 
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TABLE 5.9 

Cumulative Percent of Agricultural Floodplain 
Within Certain Absolute Error Limits for 
C-Band Simulation Results at e » 7.5® - 9.3® 


All pixels in the data base where elevation is 
> 820 feet are excluded and all non-agricultural 
pixels are excluded from analysis. 


Approx. User Accuracy Level-*- 

1 

2 

3 

4 

5 

+/- % of field capacity-*- 

50 

40 

30 

20 

10 

Moisture 

Condition 

+ 

Resolution 

-1- 






1 

f 

20 X 20 m 

96.4 

94.3 

89.9 

79.1 

52.4 


93 X 100 m 

99.6 

98,7 

96.2 

87.6 

63.8 


1 X 1 km 

100.0 

99.9 

99.4 

96.3 

76.8 

2 

20 X 20 m 

93.6 

88.5 

79.2, 

64.4 

39.6 


93 X 100 m 

98.9 

95.7 

87.6 

75.4 

47.2 


1 X 1 km 

100.0 

99.9 

98.8 

92.8 

68.7 

3 

20 X 20 m 

97.1 

95.3 

91.3 

80.9 

53.2 


93 X 100 m 

98.9 

98.0 

96.0 

89.2 

64.0 


1 X 1 km 

99.0 

98.2 

94.6 

85.9 

60.6 

4 

20 X 20 m 

97.9 

95.9 

89.8 

71.4 

39.7 


93 X 100 m 

97.2 

96.2 

93.6 

81.2 

47.4 


1 X 1 km 

89.0 

84.4 

82.7 

77.6 

50.5 
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TABLE 5.10 

Cumulative. Percent of Agricultural Floodplain 
Within Certain Absolute Error Limits for 
C-Band Simulation Results at 0 ■ U.l®-12.8® 


Mol sture 
Condition 


2C X 20 m 
93 X 60 m 
900 X 600 m 
20 X 20 m 
93 X 60 m 
900 X 600 m 
20 X 20 m 
93 X 60 m 
900 X 600 m 
20 X 20 m 
93 X 60 m 
900 X 600 m 


Approx. User Accuracy Level 

1 

2 

3 

m' 

5 

+/ 

% of Field Capacity 


50 

40 

30 

20 

10 

96.6 

94.7 

90.4 

78.3 

48.9 

98.4 

97.1 

94.7 

87.2 

58.7 

99.6 

99.2 

98.4 

94.0 

73.6 

95.0 

91.3 

82.8 

66.5 

40.3 

98.4 

96.4 

89.4 

74.6 

49.3 

99.9 

99.8 

98.6 

93.2 

73.2 

97.0 

94.5 

88.9 

75.5 

46.4 

97.7 

96.5 

92.9 

83.4 

55.4 

97.6 

96.7 

92.8 

81.9 

54.6 

96.40 

91.97 

81.1 

61.0 

34.1 

96.8 

94.5 

87.2 

70.8 

37.7 

83.9 

83.2 

82.2 

78.8 

53.7 















% of Database with Less Error 


ORIGINAL PAGE If 
OF POOR QUALITY 


100 meters 


20 meters 


Moisture Condition :01 
Angle of Incidence 7.5® to 9.3® 
183, 000 Agricultural Pixels 
on Floodplain 


Absolute Difference Between Actual and Predicted Mpc 

Figure 5.28. Soil Moisture Estimate Accuracy for Agricultural Categories on the 
Floodplain with Moititure Condition 1, oil Soil gt Field Capacity, 
and at Angles of Incidence between 7.5® and 9.3®. 






ORIGINAL PAGE IS 
OF. POOR QUALITY. 



Figure 5,29 Soil moisture estimate accuracy for agricultural categories 
on the floodplain with Moisture Condition 2, saturated soil, 
and at angles of incidence between 7.5® and 9.3®. 
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% of Database with Less Error 


ORIGINAL PAGE IS 
or POOR QUALiiy 



Figure 5.30 Soil Moisture Estimate Accuracy for Agricultural Categories on the 
Floodplain with Moisture Condition 3 , 10 Days after Thunderstorm, 
and at Angles of Incidence between 7.5^ and 9.3^. 
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% of Database with: Less Error 


ORIGINAL PAGE IS 
OF POOR QUALITY 


100 yx 

meters J 

1U 


20 meters 


Moisture Condition -.04 
Angle of Incidence 7.5° to 9.3° 
183,000 Agricultural Pixels on Floodplain 


0 20 40 60 80 100 120 140 

Absolute Difference Between Actual and Predicted Mpc 

Figure 5.31 Soil moisture estimate accuracy for agricultural categories 
on the floodplain with Moisture Condition 4, drought, and 
at angles of incidence between 7.5® and 9.3®. 
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Figure 5,32 Soil moisture estimate accuracy for agricultural categories 
on the floodplain with Moisture Condition 3, 10 days after 
thunderstorm, and with angles of incidence between 11.1° 
and 12.8®. 
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Mpc * V- 2051 plus Mp^, 785» to 965» of the agricultural floodplain Is 
accurately mapped for 0-5 cm soil moisture by the 1 km resolution 
radar, 8151 to 895* by the 100 meter resolution system, and 64X to 81 X 
by the 20-meter resolution system depending upon the moisture condition, 

The lower accuracy of the 1 km resolution in Figure 5.31 Is 
caused by the averaging of relatively high 0 ° for hard cultural targets 
over a larger area. In the case of the simulated floodplain region, 
this Is primarily related to the presence of a railroad running east to 
west and secondarily to scattered point targets such as buildings. 

The effects of resolution and moisture condition cn estimate accuracy 
are summarized In Figures 5.33 and 5.34 for the 7.5® to 9.3® angular swath 
at Accuracy Levels 3 and 4 respectively. Accuracy of moisture estimation 
using Equation 5.3 Is shown to be significantly higher on the agricultural 
floodplain than for the data base as a whole. Similar results were ob- 
tained for the 11.1® to 12.8® angular swath. 

Effects of Crop Type, Surface Roughness, and Row Direction 

The area percent of the agricultural floodplain occupied by each 
of 13 major agricultural categories of crops and surface roughness are 
given In Table 5.11. Row crops are differentiated as to row orienta- 
tion relative to radar look direction. Approximately 33^ of the area 
Is comprised of bare soil and the rest Is mainly cropped with corn, 
mi To, soybeans, and wheat. The percentage of pixels for each crop class 

A 

which have estimated moisture MpQ within +/- 20 % from the actual Mpj. 
are tabulated for each resolution and moisture condition in Tables 5.12 
and 5.13. For the 20-meter resolution (and to a lesser extent the 100 meter 
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Resolution (Meters) Resolution (Meters) 

Figure 5.34 Percent of data base where estimated moisture Mpc is within ±2056 of 
actual moisture Mpc. Angle of incidence e is 7.5° to 9.3° relative 
to nadir. Accuracy level 4 is approximately ±20° of field capacity. 


TABLE 541 


Am Percent of the Agricultural Floodplain for 
Crop Categories; Floodplain Is Defined by 
Elevation s 820 Feet 


Crop Category 

Area 

Percent 

Bare Soil - Smooth 

15.20 

Medium Rough 

10.26 

Rough 

7.09 

All Roughnesses 

32.55 

Pasture 

5.13 

Mown Pasture 

0.51 

Alfalfa 

5.54 

Soybeans - Parallel 

5.54 

Perpendicular 

5.36 

All Soybeans 

10.90 

Milo - Parallel 

6.96 

Perpendicular 

7.50 

All M1lo 

14.46 

Corn - Parallel 

9.56 

Perpendicular i 

11.88 

All Com 

21.44 

Wheat 

9.47 

All Crops and Pasture 

67.45 
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TABU 5.12 


-f/- 20 % Accuv*acy Levels for Crop 
Categories on the Floodplain for the 
7.5® to 9,3® Angular Swath 


Percent of Category Within f/- 20^ Estimate Accuracy 


Moisture Condition 

1 

2 

3 

4 

Resolution 














Crop 4^ 

% of 
Flood- 
plain 4 

20m 100m 1km 

20m loom 1km 

20m 100m Ikm 

20m 100m 1km 

All Classes 

100 .o' 

79 

88 

96 

64 

75 

93 

81 

89 

86 

71 

81 

78 

Bare * Smooth 

15.2 

31 

91 

97 

82 

92 

97 

72 

86 

81 

56 

77 

83 

Bare - Medium 
Rough 

10.3 

86 

95 

95 

85 

95 

96 

82 

93 

90 

80 

93 

85 

Bare Rough 

7.1 

86 

93 

94 

85 

91 

92 

87 

94 

82 

87 

95 

85 

Pasture 

5.1 

68 

76 

97 

50 

67 

91 

69 

77 

87 

68 

79 

71 

Alfalfa 

5.5 

79 

82 

86 

54 

71 

98 

81 

80 

71 

76 

75 

49 

Soybeans // 

5.5 

88 

96 

97 

84 

94 

89 

90 

95 

84 

90 

93 

72 

Soybeans X 

5.4 

90 

95 

97 

86 

94 

98 

87 

92 

89 

87 

89 

80 

Milo // 

7.0 

45 

55 

93 

5 

12 

71 

78 

88 

84 

84 

88 

64 

Milo i 

7.5 

57 

68 

98 

9 

17 

91 

82 

91 

91 

76 

83 

74 

Corn // 

9.6 

88 

96 

98 

72 

83 

90 

86 

93 

93 

61 

72 

83 

Corn 1 

11.9 

87 

95 

99 

75 

86 

94 

79 

87 

87 

50 

57 

81 

Wheat 

9.5 

85 

93 

99 

57 

74 

99 

86 

93 

89 

81 

91 

79 
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TABLE 5.13 

± 20% Accuracy Levels for Crop Categories on the Floodplain 
for the 11.1® to 12.8° Angular Swath 


I 


Moisture Condition 


Resolution 

s 

iOiD<Vir-0>^OU>hoQrs. 

100 !■ 

71 

25 

87 
95 
84 
79 
95 

88 
87 
67 
74 
48 
92 

c 

o 

CVJ 

61 

14 

74 

86 

76 

77 
89 
79 
74 
55 
62 
43 
79 

n 

1 Km 

82 

78 

87 

81 

77 

67 

80 

87 
81 

88 
91 
82 
80 

loom 

83 

58 

93 

94 
74 
77 
94 

91 
90 

92 

93 
81 
90 

20m 

76 

so 

82 

87 
67 

77 

88 
82 
82 
84 
84 
70 
82 

CvJ 

1 Km 

93 

96 

98 

98 
83 
91 
93 

99 
77 
93 
91 

96 

97 

E 

s 

r* 

75 

91 

62 

7S 

67 

79 

94 

93 

18 

28 

87 

90 

89 

20m 

67 

82 

48 

56 

60 

75 

87 

85 

16 

32 

77 

82 

82 

- 

1 Km 

94 
93 

95 

95 
89 
79 
93 
99 
93 
98 

96 
98 
96 

100m 

87 

84 

91 

94 
68 
75 

96 
93 
73 

97 

95 
93 
89 

20m 

— 

78 

74 

77 

80 

64 

76 
88 
85 
64 

77 
87 
83 
82 

Percent 

of 

Floodplain 

100.0 

15.2 

10.3 

7.1 

5.1 
5.5 
5.5 

5.4 
7.0 

7.5 

9.6 
11.9 

9.5 

Crops 

All classes 

Smooth Bare 

Medium Rough 
Bare 

Rough Bare 
Pasture 
Alfalfa 
Soybeans // 
Soybeans X 
Milo // 

Milo 1 
Com // 

Com X 
Wheat 
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resolution) most of the estimate error Is concentrated within 
the categories of mllo, wheat, alfalfa, and com for saturated 
moisture conditions and within the categories of smooth bare soil and 
corn for very dry moisture conditions; and the accuracy of predic- 
tion Is observed to decrease at higher incidence angles (Table 5.13). 

For the 20-meter resolution, the concentration of a large percent- 
age of net estimate error within a relatively few crop classes is 
to be expected from the use of a general estimation algorithm 
such as Equation 5.3 unless that algorithm Is weighted by the ex- 
pected percentage of crop classes for a given region. Importantly, 
the need for a regional and/or time-dependent general estimation al- 
gorithm appears to be mitigated in large part by the natural averaging 
of local canopy effects when using a coarser resolution system. Thus, 
the inter-class variance in estimate error is seen to be very low for 

the 1 km resolution at all simulated moisture conditions in Tables 
5.12 and 5.13. The only exception to the preceding statement is 
alfalfa, where the large estimate error for Moisture Condition 4 from 
the 1 km resolution is an artifact of its close proximity to the 
previously mentioned railroad. 

Because of the high moisture estimate accuracy attained for specific 
agricultural classes by the foregoing approach using Equation 5.3 as the 
estimation algorithm, the division of the data base into bare soil and 
vegetation-covered classes for differential estimation of soil moisture 
using Equations 5.4 and 5.5 was deemed unnecessary. For the simulated 
conditions, it seems unlikely that the use of a prioH information other 
than the location of water bodies and cultural features would signifi- 
cantly improve estimate accuracy and be worth the added complexity in 
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processing. The distribution of prediction error between target classes 
observed in Tables 5,12 and 5.13 indicates that surface roughness and 
specific crop canopy type are nwre significant determinants of net 
prediction accuracy than a dichotomous division of the agricultural 
terrain into bare and vegetated classes, 

6.0 CONCLUDING REMARKS 

The validity of results derived from any simulation of a real-world 
situation Is inherently limited by the degree of realism incorporated in 
the assumptions and models used In generating the simulations. In the 
present study, the aspect of realism is governed by three types of fac- 
tors: (a) geometrical factors associated with the image formation pro- 
cess of a sidelooking imaging radar, (b) signal fluctuations due to 
Rayleigh fading (speckle effect) in radar images that contain a small 
number of independent samples per pixel element, and (c) the models 
characterizing the backscatter behavior of the various target classes 
contained in the simulated scene. The first two types of factors are 
well understood/, therefore, it was possible to incorporate them in the 
simulation algorithm v^/ith a high degree of accuracy. This accuracy 
was verified in a separate study through comparisons of simulated SAR 
images with actual images of the same scene. The degree of realism 
associated with the third type of factor— the backscatter models— is 
very good, if considered in statistical terms. That is, the functions 
used to describe the dependence of the backscattering coefficient a® 
on the local angle of incidence and moisture content (for agri- 
cultural categories) are based on statistical regressions applied to 
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experimental data obtained over the past eight years. In the majority of 
cases, the correlation coefficients associated with these regressions are 
greater than 0.8, which means that the regressions account for the major- 
ity of the observed variation In but some variation, part of which Is 
due to measurement error, remains unaccounted for. Considering that the 
overall study Is, to some extent, statistical In nature, the approach used 
above for modeling the backscatter behavior Is certainly justified. 

The major conclusions derived from this study are; 

(1) Among the three resolutions considered (20 m, 100 km, and 1 km), 
soil moisture content was estimated with the highest accuracy by the 1-km 
resolution radar images for the relatively wet soil conditions (Cases 1 
and 2 of Table 4.7), comparable levels of accuracy were provided by the 
100 m and 1 km resolutions for Case 3 (10 days after the thunderstorm), 
and the 100 m resolution provided significantly better results than the 

1 km resolution for the drought case (30 days after the thunderstorm). 

(2) Based on the above results. It appears that a spatial resolu- 
tion between 100 m and 1 km would provide optimum performance over the 
various soil moisture conditions. Narrowing this range down to a sped- 
fic value is the object of a further phase of this investigation. 

(3) This study was performed for a test site in eastern Kansas, 
where dry-land farming practices prevail. For such regions, the peri- 
odic structure associated with soil surfaces of row crops exercises a 
minor influence on the scattering behavior as a function of radar look- 
direction relative to row-direction (for the sensor parameters specified 
In this study). In regions where irrigation practices are common, the 
row structures usually have larger amplitudes, which would lead to greater 
ambiguity In the soil moisture estimation process for like-polarized 
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radar configurations.. Research conducted to date Indicates that the | 

cross-polarizod scattering coefficient Is significantly less sensitive 
to row direction and, therefore, should be preferred for napping soll- 
molsture content In Irrigated regions. The major drawback to the use 
of cross-polarization Is the fact that more transmitter power Is 
required than for the like-polarized case. However, since the needed 
resolution Is of the order of hundreds of meters. It may be possible 

I 

to configure a cross-polarized space radar system with | 

existing technology. | 


I 
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